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Abstract 

We study the infinite horizon Linear-Quadratic problem and the asso- 
ciated algebraic Riccati equations for systems with unbounded control ac- 
tions. The operator-theoretic context is motivated by composite systems 
of Partial Differential Equations (PDE) with boundary or point control. 
Specific focus is placed on systems of coupled hyperbolic/parabolic PDE 
with an overall 'predominant' hyperbolic character, such as, e.g., some 
models for thermoelastic or fluid-structure interactions. While unbounded 
control actions lead to Riccati equations with unbounded (operator) co- 
efficients, unlike the parabolic case solvability of these equations becomes 
a major issue, owing to the lack of sufficient regularity of the solutions to 
the composite dynamics. In the present case, even the more general the- 
ory appealing to estimates of the singularity displayed by the kernel which 
occurs in the integral representation of the solution to the control system 
fails. A novel framework which embodies possible hyperbolic components 
of the dynamics has been introduced by the authors in 2005, and a full 
theory of the LQ-problem on a finite time horizon has been developed. 
The present paper provides the infinite time horizon theory, culminating 
in well-posedness of the corresponding (algebraic) Riccati equations. New 
technical challenges are encountered and new tools are needed, especially 
in order to pinpoint the differentiability of the optimal solution. The 
theory is illustrated by means of a boundary control problem arising in 
thermoelasticity. 



Introduction 

Historical background and motivation. The theory of the optimal con- 
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trol problem with coercive quadratic functionals for abstract linear systems 
y' = Ay + Bu — in spaces of finite or infinite dimensions — is intrinsically linked 
to well-posedness of the corresponding algebraic/differential Riccati equations. 
As it is well known, by solving these celebrated equations one obtains the op- 
erator (a matrix, in the finite dimensional case) which occurs in the feedback 
representation of the optimal control as well as in the quadratic form which 
yields the optimal cost. 

In an infinite dimensional setting, if B is a bounded control operator from 
the control space U to the state space Y, then a complete theory of the Riccati 
equations has been developed, where the analogs of all the key properties known 
in the classical finite dimensional context hold true. Namely, the unique solu- 
tion to the Riccati equation provides the optimal cost operator, the feedback 
representation of the optimal control holds true, and the closed loop equation 
yields the synthesis of the optimal control; see [TU] for pioneering results in this 
direction, [TT]. [T3] and [T5]. 

When it comes to the abstract representation of initial/boundary value prob- 
lems for evolutionary Partial Differential Equations (PDE) in a bounded domain 
f2 C M", with control action exercised on the boundary dfl or pointwise in the 
interior of the domain, the major technical challenges come from the unbound- 
edness of the control operator B which naturally arises from the modeling of 
the PDE problem. This, in turn, results in the possible unboundedness of the 
gain operator B*P which occurs in the quadratic term of the Riccati equations. 
In fact, this operator may be even not densely defined. Indeed, this 'pathology' 
is shown to happen even in the case of simple hyperbolic equations with point 
control; see [37], [35] . 

In this respect, the regularity properties of the pair (A, B) which describes 
the free dynamics and the control action or, more specifically, of e At B (where 
e At is the semigroup governing the free dynamics) are absolutely central. And 
in fact, historically speaking, the first extensions of the LQ-problem theory 
for systems with bounded control operator to the present setting pertained 
to parabolic-like dynamics, where the underlying semigroup e At is analytic. 
(Canonical illustrations are the heat equation with boundary or point control, 
as well as structurally damped plate equations, and certain thermoelastic sys- 
tems.) Analyticity made it possible the full development of a theory of the 
Riccati equations, with unbounded gains as well as without additional regu- 
larity properties of the optimal cost operator. These results date back to the 
beginning of the eighties ([25]); see [26] for a concise treatise, or, alternatively, 
the extended monographs [13] and [30] . The analytic theory of Riccati equations 
with non- autonomous coefficients has been developed in [3J, [3], [5J|S]. 

In contrast, as it was pointed out already, the lower regularity exhibited by 
e At B in the case of hyperbolic- like equations with boundary /point control is 
in general insufficient to guarantee the sought-after properties recorded above. 
(Canonical examples are provided by PDE which display some kind of 'wave 
propagation', even in one dimension, under unbounded control actions.) Thus, 
a requirement of an appropriate smoothing effect of the observation operator is 
called for, in order to establish well-posedness of the differential Riccati equa- 
tions. See, again, [55], [T3], [3D], and the more recent [T2"] . 

The first break with analyticity and additional assumptions on the observa- 
tion operator was introduced by G. Avalos in 1996 in the study of PDE models 
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that arise in structural acoustics; see |S]. In these problems the equation 
for the acoustic waves is coupled on an interface with a parabolic-like equa- 
tion describing the structural vibrations of an elastic wall and therefore the 
systems under consideration comprise both hyperbolic and parabolic dynamics. 
The study carried out in [8] pinpoints and establishes a fundamental (singular) 
estimate for the norm of the operator e At B in the space C(U, Y); this is rec- 
ognized as the actual mathematical property which by itself — in the absence of 
analyticity — is sufficient to ensure boundedness of the gain operator. It is worth 
noting that while the aforesaid 'singular estimates' are intrinsic to control sys- 
tems whose free dynamics is governed by an analytic semigroup e At , in the case 
of general Co-semigroups these estimates correspond to suitable (non trivial) 
interior regularity results for the solutions of the uncontrolled system of PDE. 

This path of investigation has been followed later on, bringing about a the- 
ory of the LQ-problem on both a finite and infinite time horizon, while other 
significant PDE systems yielding 'singular estimates' have been discovered and 
analyzed; see [23j [24] , [23 [29], and [31] (dealing with the Bolza problem). It 
turned out that the class of control systems characterized by this property — 
although inspired by the model for acoustic-structure interactions studied in 
[5] — covers a variety of systems of coupled hyperbolic/parabolic PDE, includ- 
ing different structural acoustic models, thermoelastic plate models, composite 
(sandwich) beams models, and others; see, e.g., [23], [29] Section 4], [24] Lec- 
ture III, Part II], [TB], [32] ■ (See also [T3] for further applications of singular 
estimates to the study of semilinear evolution equations with nonlinear bound- 
ary terms.) 

In fact, interactions between distinct physical phenomena widely occur in 
both nature and technology. To name a few relevant ones, we just recall — 
besides thermoelastic and acoustic-structure interactions — fluid-solid or magne- 
toelectric interactions. The modeling of such interactions leads in a natural way 
to composite systems of PDE comprising components which display different 
dynamical behaviours. In a mathematical setting, such evolutionary systems 
of coupled PDE can be represented in a simple form by the following abstract 
system, which also takes into account the possible influence of control actions: 

y t = Axy + B lU + Ci(y,z) y\ t=Q = y° G H x 
Zt = A 2 z + B 2 v + C 2 (z,y) z\ t=0 = z° e H 2 ' 

where H± and H 2 are appropriate Hilbert spaces, A\ and A 2 are infinitesimal 
generators of Co-semigroups e Ait on Hi, i = 1,2, respectively, with e.g. e Alt 
assumed to be analytic. The control operators Bi, i = 1,2, are unbounded 
operators acting from the control spaces f/j into Hi (or from Ui into [J9(yl*)]'); 
the coupling between the two different dynamics may occur as well by means of 
unbounded operators C. 

It is then natural to seek the proper set of assumptions on the overall dy- 
namics operator in order to fully exploit the diverse features of the components 
of the system. The ultimate goal is to produce a rigorous theory which yields 
the feedback synthesis of the optimal control by means of a solution to the (non- 
linear) operator Riccati equation, which must be shown to be well posed. The 
unboundedness of the control operator, while being a natural and prominent 
feature of the problem, provides the main mathematical challenge to achieve a 
proper definition of the gain operator, which may not have a sufficiently rich 
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(and hence, meaningful) domain. On the other hand, the framework of control 
systems characterized by singular estimates turned out to be too rigid and fail- 
ing to encompass many 'mixed' dynamics that are of fundamental importance 
in the applications. 

A prime illustration of this fact is a benchmark thermoelastic plate model 
studied in [16 . The distinct control-theoretic properties exhibited by its ab- 
stract representation, according to different sets of mechanical/thermal bound- 
ary conditions, revealed that the singular estimate which pertains to its parabolic 
component (that is the equation for the temperature distribution in the plate) 
does not always fully 'propagate' to the entire system. Under certain boundary 
conditions, the coupling brings about a more involved behaviour of the (con- 
trolled) coupled dynamics. More precisely, in contrast with the case of hinged 
boundary conditions, when the system is supplemented with clamped (mechani- 
cal) boundary conditions and is subject to Dirichlct thermal control, the coupled 
dynamics precludes the validity of a singular estimate (for the norm of the oper- 
ator e At B, as a bounded operator from the control space U into the state space 
Y); see [S] and [J. 

It was indeed this initial/boundary value problem for the system of ther- 
moelasticity to provide our original motivation for introducing in [2] a novel 
class of abstract linear control systems — which includes the one discussed above 
as a very special case. This class is characterized by a suitable decomposition 
of the operator e At B = F(t) + G(t), where the only component F(t) satisfies 
a singular estimate, while the component G(t) exhibits appropriate regularity 
properties which account for the (predominant) hyperbolic character and the 
parabolic one; see Hypothesis 2.2] for full details. 

The goal was then to develop a sufficiently general abstract set-up which could 
encompass significant interconnected PDE systems comprising both hyperbolic 
and parabolic components, and yet lacking the overall 'parabolic-like' character 
disclosed by a singular estimate. Moreover, we aimed at developing a corre- 
sponding theory of the quadratic optimal control problem which escaped re- 
strictions on the observation operator in order to achieve well-posedness of the 
Riccati equations, unlike the case of purely hyperbolic problems. 

Thus, a theory of the LQ-problem for the class of systems described above 
(and mathematically defined by Hypothesis 2.2 in [2]) in the finite time horizon 
case has been developed in [2J. The major novel features of this optimal control 
theory are shortly outlined below; the reader is referred to [H Theorem 2.2] for 
a complete description of the obtained results. 

First of all, in contrast with previous theories the gain operator B*P{t) is 
only densely defined in the state space Y (it is bounded on T>{A e ), with arbitrarily 
small e) . This turns out to be sufficient to obtain well-posed Differential Riccati 
equations 

^(P(t)x, z) Y + (A*P(t)x, z) Y + (P(t)Ax, z))Y - (B*P(t)x, B*P{t)z) v 
at 

+ (R*Rx, z) Y = for any x, z G V{A). 

Secondly, solvability of the corresponding quadratic optimal control problems 
follows without assuming smoothing properties of the observation operator R. 
More precisely, the weak requirement on R — that is the same as condition (jl.8|) 
of Hypothesis 11.41 — allows bounded operators which, roughly, just 'maintain 
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regularity', such as the identity operator. Consequently, natural quadratic func- 
tionals such as the ones which involve the integral of the physical energy of the 
system are allowed for the first time. 

Model-specific analyses needs to be carried out on diverse PDE problems, 
in order to establish the regularity of boundary traces that is required in order 
to apply the introduced abstract theory (of the LQ-problem on a finite time 
interval). This has been achieved, indeed, in the case of boundary control 
problems for acoustic-structure and fluid-structure interactions as well; see [15] 
and [H1[TB], respectively. 

Present work. The purpose of this paper is to complement the LQ-theory 
described above by developing a complete infinite time horizon analysis. The 
task is not straightforward owing to a natural mixing of singularities occurring 
in short and long time. The interplay between the long time stability for the 
forward problem and the short time development of singularities for the adjoint 
problem lie at the heart of the problem. 

We emphasize at the outset that we establish solvability of the optimal con- 
trol problem, as well as well-posedness of the corresponding Algebraic Riccati 
equations, under minimal assumptions on the operators involved. More pre- 
cisely, we work on the very same abstract framework set forth in the finite 
time horizon case, that is the one defined by [2j Hypothesis 2.2] for the oper- 
ator e At B and the observation operator R which occurs in the cost functional 
(|f .3[) . The finite time horizon scenario is only complemented with two natural 
stability requirements pertaining to the underlying semigroup e At and the com- 
ponent F(t) arising from the decomposition of e At B; see Hypothesis 11.11 and 
Hypothesis 11.4( 1) in the next Section. 

As a general strategy for the proof, we start with the variational approach 
which has been pursued in the work of the authors of [3U], since their ear- 
lier studies of parabolic PDEs with boundary control; see, e.g., [25] ■ Accord- 
ingly, the program evolves through the following steps, (a) The existence of a 
unique optimal pair (u, y) follows by convex optimization arguments, (b) The 
Lagrange-dual multipliers method yields the optimality condition for the opti- 
mal pair. Then, the optimal pair (u, y) is characterized in terms of the data of 
the problem, (c) An operator P is introduced, defined in terms of the optimal 
state y(t; yo) = $(i)yo; this operator is ultimately shown to satisfy the Algebraic 
Riccati Equation (ARE). 

However, although the general philosophy of the aforesaid approach applies 
as well to the present abstract class of boundary control systems, the regularity 
specific to the operator e At B — equivalently, to the corresponding input-to-state 
map L — brings about novel technical challenges in order to establish that the 
ARE is well-posed. Delicate issues are encountered when we are 

(i) to give a meaning to the gain operator B*P, for which boundedness on 
the state space does not hold: hence, it is sought (and in fact established) 
on a dense subset; even more, when we are 

(ii) to pinpoint the regularity of the map t <— > $(t)B, which plays a central role 
in the study of the differential properties of the optimal state semigroup 
&(t) — that is eventually shown to be differentiable for t > on T>{A). 

We note that T>(A) is not a natural domain of the strongly perturbed evolution 
<£>(£) and therefore the aforementioned differentiability is far from expected. 
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This differential property is a consequence of an appropriate (time and space) 
regularity result obtained for the map t i-> &(t)B, which in turn hinges on the a 
suitably developed operator perturbation theory applied to the original operator 



It is important to emphasize that while the former issue (i) has its counter- 
part in the finite time horizon theory for the present class of dynamics (see [5] ) , 
(ii) requires novel developments not only with respect to previous theories of 
the LQ-problem (namely, for different classes of boundary control systems, such 
as the one characterized by singular estimates), but also with respect to the 
finite time horizon case for the present class of dynamics. For our developments 
we introduce and employ an appropriate class of weighted function spaces (see 
(I1.15|) ), whose central role will become fully apparent in Section l4~2l 

Paper outline. We conclude the Introduction with a brief overview of 
the paper. In Section[T]we introduce the class of linear control systems, charac- 
terized by Assumptions 11.1111.41 for which we develop the infinite time horizon 
optimal control theory. The major statements of this theory are collected in 
Theorem 11.51 Few remarks about the notation are given at the end of the 
Section. 

In Section [2] we derive a complex of regularity results which concern the 
operator B*e A *, as well as the components F and G of its decomposition. 
In addition, we develop a full regularity theory for the input-to-state map L 
defined by (12.101) and its adjoint L* (Proposition 12. 6[ Proposition 12.71) . These 
results constitute the fundamental basis for the developments of the subsequent 
sections. 

In Section [3] we introduce the optimal cost operator P and show that the 
gain operator B*P is bounded on T>(A e ) (Theorem I3.5|) . The validity of the 
feedback representation of the optimal control is established here. 

In SectionS]we prove well-posedness of the algebraic Riccati equations. The 
technical basis for the corresponding Theorem 14. 161 is found in Proposition ^. 131 
and its Corollary 14.141 which in turn follow from the novel result of Theo- 
rem [47T2J The preliminary description of the domain V(Ap) of the optimal 
state semigroup, along with the detailed information provided by Proposition l4.5l 
(whose proof relies on the theory of interpolation spaces) , will also prove essen- 
tial in showing Corollary 14. 141 

In Section [S] we briefly illustrate the applicability of the obtained infinite- 
time horizon theory through the analysis of a natural optimal control problem 
for the thermoelastic system which motivated and initiated our former theory 



1 The mathematical problem. Main results 

Let Y and U be two separable Hilbert spaces, the state and control space, 
respectively. We consider, on the extrapolation space [D(A*)}' , the abstract 
control system 



e At B :U-¥ [D(A*)}' . 



of 0. 




t > 



(1.1) 



under the following basic Assumptions. 
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Assumption 1.1. Let Y , U separable complex Hilbert spaces. 

• The closed linear operator A : T)(A) C Y — > Y is the infinitesimal gen- 
erator of a strongly continuous semigroup {e At }t>o on Y , which is expo- 
nentially stable; namely, there exist constants M > 1 and tu > such 
that 

\\e At \\ C (Y) < Me~ ut Vi>0; (1.2) 

then A' 1 G £(Y); 

• B : U — >• is a bounded linear operator; equivalently, A~ 1 B G 
C(U,Y). 

To the state equation (jl.ip we associate the quadratic functional 

/>oo 

J(u)= / (\\Ry(t)\\% + \Ht)\\u)dt, (1.3) 
Jo 

where Z is a third separable Hilbert space — the so called observation space 
(possibly Z = Y) — and at the outset the observation operator R simply satisfies 

R€C(Y,Z). (1.4) 

Remarks 1.2. (i) Since by Assumption 11.11 the semigroup e At is uniformly 
stable, then —A is a positive operator and the fractional powers (— A) a , a € 
(0, 1), are well defined. In order to make the notation lighter, we shall write A a 
instead of (— A) a throughout the paper. 

(ii) We note that the functional (ll.3[) makes sense at least for u = 0. This again 
in view of the exponential stability of the semigroup e At (Assumption II. ip . 
which combined with (|1.4[) ensures Ry(-,yo;0) G L 2 (Q,oo;Y). 

(iii) The analysis carried out in the present paper easily extends to more general 
quadratic functionals, like 

/>oo 

J(u)= / (\\Ry(t)\\%+\\Ru(t)fu)dt, 
Jo 

provided R is a coercive operator in U . We take R = I just for the sake of 
simplicity and yet without loss of generality. 

The optimal control problem under study is formulated in the usual way. 

Problem 1.3 (The optimal control problem). Given yo 6 Y, seek a control 
function u G L 2 (0, T; U) which minimizes the cost functional (|1.3[) . where y(-) = 
y(- ;yo,u) is the solution to corresponding to the control function u(-). 

Aimed to pursue the study of the infinite horizon problem for the abstract 
class of boundary control systems first introduced in [2] , we assume throughout 
the paper that the dynamics, control and observation operators are subject to 
the following conditions. 

Assumptions 1.4. The operator B*e A 1 can be decomposed as 

B*e A " t x = F(t)x + G(t)x, t>0, xeV(A*), (1.5) 

where F(t) : Y — > U and G(t) : T>{A*) —±U,t> 0, are bounded linear operators 
satisfying the following assumptions: 
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(i) there exist constants 7 € (0, 1) and rj, N > such that 

\\F(t)\\ c{ Y U) <Nt-'e-^ Vi>0; (1.6) 

(ii) there is a time T > such that the operator G(-) belongs to C(Y, L p (0, T; U)) 

for all p € [1, oof; 

(iii) there exist T > and e > such that: 

(a) the operator G{-)A*~ t belongs to C(Y, C{[0, T]; U)), and in particular 

K T := sup \\G{t)A*- f -\\ c{YiU) <oo; (1.7) 
te[o,T] 

(b) the operator R*R belongs to C{V(A e ),V(A* e )), i.e. 

\\A* t R*RA- £ \\ c{Y) <c<oo; (1.8) 

(c) there exists q G (1,2) smc/i that the operator B*e A ' A* t , which is de- 
fined in T>{ A* has an extension which belongs to C(Y, L q (0,T; U)) 

1.1 Statement of the main results 

The main result of this paper is the following Theorem, which collects the most 
significant statements which pertain to the solution of the optimal control prob- 
lem, as well as to well-posedness of the corresponding algebraic Riccati equa- 
tions. 

Theorem 1.5. Consider the optimal control Problem \1.3[ under the set of 
Huvotheses ll. lf[1.4\ Then, the following statements are valid. 

51. For any yo £ Y there exists a unique optimal pair (&(•),$(•)) for Prob- 
lem \1.3l which satisfies the following regularity properties 

ue (J L p (0,co;£7), y G C 6 ([0, 00); Y) H [ \J L p (0,oo ; y)]. 

2 < p < oo 2<p< 00 

(1.9) 

52. The operator $(<), t > 0, defined by 

$(t)y :=y(t)=y{t,y ;u) (1.10) 

is a strongly continuous semigroup on Y , t > 0, which is exponentially 
stable; namely, 

3Mi>l,wi>0: < Mi e"" 1 ' Vi > . (1.11) 

53. The operator P G C(Y) defined by 

/>oo 

Px := / e A ' t R*R<5>{t)xdt xeY (1.12) 
J 

is the optimal cost operator: namely, 

/•OO 

(Px,x) Y = / (\\Ry(fr,x))\\ 2 z + \\u(t;x)\\l)dt, e y , 

Jo 

which also shows that P is (self-adjoint and) non-negative. 
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S4- The gain operator B*P belongs to C(T>(A e ),U) ; namely, it is just densely 
defined on Y and yet it is bounded on T>(A e ). 

55. The infinitesimal generator Ap of the (optimal state) semigroup $(t) co- 
incides with the operator A(I — A~ 1 BB*P), on the domain 

V{A P ) C {x e Y : x - A- 1 BB*Px E V{A)} 

i rt 

C {x e Y : 3w-lim f ^ 0+ - / e A ^~ T) A' 1 BB* P<S>(t)x dr 
t Jo 

56. The operator e At B, defined in U and a priori with values in [D(A*)}' , is 
such that 

e s e A B e C(U, L»(Q, oo; [V(A* e )]')) Vp € [1, l/ 7 ) (1.13) 

for all 5 € [0,wA r/); almost the very same regularity is inherited by the 
operator $(t)_B; 

e s <f>(-)Be £{U, L p (0, 00; [V(A* £ )}')) Vp g [1, 1/7) , (1.14) 
provided S S [0, cj A 77) is sufficiently small. 

57. The optimal cost operator P defined in (|1.12[) satisfies the following addi- 
tional regularity properties: 

P 6 C(D(A P ),V(A*)) n C(V(A), V{A P )) ; 

moreover, P is a solution to the Algebraic Riccati equation corresponding 
to Problem ] 1.31 that is 

(Px,Az) Y + (Ax,Pz) Y - (B*Px,B*Pz)u + {Rx,Rz) z = 

for any x, z G T>(A), 

to be interpreted as 

{A*Px, z) Y + {x, A*Pz) Y - (B*Px, B*Pz) v + (Rx, Rz) z = 

when x, z € V(Ap). 

58. If x 6 T>{A e ), the regularity (|1.9p of the optimal pair is improved as follows: 

y 6C b ([0,oo);P(A e ))n [ |J L p (0, 00; £>(A e ))] , 

1<;j<oo 

fi e C 6 ([0, 00); 17) n [ (J L p (0,oo;f/)] . 

l<p<oo 

TTie following (pointwise in time) feedback representation of the optimal 
control is valid for any initial state x 6 Y : 

u{t,x) = —B*Py(t,x) for a.e. t £ (0, 00). 
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1.2 Notation 

Inner products in Hilbert spaces X (Y and U, in practice) will be denoted by 
(■, -)x throughout the paper. The subscripts will be omitted when no confusion 
arises. Instead, the symbol {-, -}y will denote a duality pairing of V with V; 
V = V(A* e ) will occur most often. 

We shall utilize L p -spaces with weights, defined (in the usual way) as follows: 



where g : (0, oo) — > M is a given (weight) function. We will more specifically 
utilize exponential weights such as g(t) — e St ; to simplify the notation we will 
write 



2 The input-to-state map: relevant regularity 
results 

We begin by providing a series of regularity results, concerning first the operator 
B*e A * (or one of its components) and then the input-to-state map L defined 
by (|2.10[) (and its adjoint L*). These results constitute the first consequences 
of the abstract Assumptions 11.41 as well as the fundamental basis for the more 
challenging developments of the subsequent sections. 

2.1 Preliminary results 

Proposition 2.1. For each 5 € [0, ui A 77 [ and p £ [1, 00 [ the map t 1— > e 5t G{t) 
belongs to C(Y, LP (0,oo; U)). 

Proof. Let T > be such that Hypotheses ll.4f ii) holds. By hypothesis ll.4f i). 



L|(0, 00; X) := {/ : (0, 00) -> X , g(-) /(•) 6 D>(0, 00; X)} , 



LftO, 00; X) := {/ : (0, 00) X , e 5 ' /(•) £ ^(0, 00; X)} 



(1.15) 



we have 



F(-) G£(F,L p (e,T/2;[/)) Ve g]0, -] . 



Hence 




T/2 



(2.1) 



Now we can write: 





(fc+l)T/2 



e* pt ||B*e j4 **x||P di 



k=1 JkT/2 



OO nj 1 




E 

fe=i 
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using (|2.1[) . we deduce for each x £ 2?(A*) 

,.00 oo 

/ ||e* t B*e A * t o;||P dt < ce^ T V e W*-i>5 U^'^i < 

OO 

< ce* T M" ^V^- 5 ^- 1 )^^ = c(p J T)||a;[|?. . 
fe=i 



By density, this shows that the map 1 1-> e <5t i?*e j4 * has a continuous extension 
from Y to L p (T/2, oo; £7). On the other hand, Hypothesis II. 4f i) implies that 

r OO / r, \ IP poo 

/ ||e 5t F(i)*||£ <tt < W - / e-M^I^I* dt < c(p, T) H^f, , 



'T/2 \ J / JT/2 

so that 



/■OO 

/ lie 54 


B *e A " t -F(t) 


IT/2 





We^G^xWldt 

T/2 

As, by Hypothesis [L^ii) , 

T/2 || e 4 *G(f)a:||* di < e 5 ^ / T/2 \\G{t)x\%dt < c(p, T)\\x\\ p Y , 



we conclude that t 1-> e St G(t) is continuous from Y to £ p (0, oo; {/). □ 

Proposition 2.2. for eac/i 5 £ [0, lu A rj[ the map t h-> e <5 *_B*e A **v4* € /ias an 
extension which belongs to C(Y, L q (0, oo; £/))• 

Proof. Let T > such that Hypothesis II .4f iii) (c) holds. We can write, for each 
x e V(A* e ), 

poo oo p(k+l)T 

/ \\e St B*e A ' t A* £ x\\ q u dt = V / e 5 «* ||B*e A **A* e x|| dt = 

Jo fc=0"'* r 

oo „x 

= J- / e ff «t'+* T 5||fl*e A *M* 6 e A, * T a:|^d«. 
As e A ' kT x G £>(A* C ) whenever a; £ 2?(A* e ), we get by Hypothesis OJiii) (c) 

r oo oo .T 

/ ||e 4t B*e A *'A* e a;||« ds < e^ T £ e^ feT / ||B* e A * s A* e [e A * fcT a;] ||« ds < 



< ce 5qT Y, Mq ^' S)9kT \\ x]]S - 



iy ■ 

fe=0 



□ 



Proposition 2.3. For each 6 £ [0,o> A 77 [ and p £ [l,l/7[ the map t H» 
e St G{t)f{t) belongs to L p (0,oo;U) whenever f £ 1/(0, oo; 2?(A* e )), witfi 

— - — <r<oo; (2.2) 
1 - TP 

m addition, 

||e' 5 'G(-)/(-)IUp(0,oo:C/) < Cp||/||ir(0 lOO;2 3(A» = )) • (2.3) 

T7ie constraints assumed on the exponents p and r are sharp. 
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Proof. Let / 6 L r (0, oo;T>(A* t )) be given, initially with 1 < r < oo. We take a 
representative of f(t), defined for almost any t £ [0, T], and derive a preliminary 
estimate of ||e*-G(-)/(-)ll in £ P (0,T;[/), that is 



\\e st 



G{-)f{-Wudt<K p T / e^||A*7WII?-^ 
Jo 

<K^j\ 5t ^-^dtf r ' P)/T 

= ^11/11^(0^.^(^.6)) 



li r (0,oo;-D(A*«)) 



(2.4) 



where we utilized first Hypothesis ll.4t iii) (a) and then we applied the Holder 
inequality with exponents r/(r — p) and r/p, taking r > p > 1. 

On the other hand, owing to Hypothesis 11.4( 1). we immediately find the 
estimate 

T pT -( V -S)pt 

e^\\F(t)f(t)\\ P udt< N p J Q — \\f(t)\\ p Y . (2.5) 

Thus, in order to render the integrand in the right hand side of (12.51) summable, 
we take p G [1, 1/7) first, and apply once more the Holder inequality, this time 
with exponents s/p and s/(s — p), with s = pr/(r — p), thus obtaining 



T 



-(»7— ff)at sp/s 

e°n\F(t)f(t)Kdt<NPl I " .... df 



ip 



tls imizrovx^A..)) (26) 

= ^pll/llz,r(o,oo;X>(A» e )) ' 

We note that in (|2.6p we used s/p > 1, along with 57 = prj/(r —p) < 1, which 
readily yields the lower bound in (|2.2I) . Combining (|2.6p with (|2.4p we find 

e «p*|| fl *eA-* /(t) ||P dt^jjp f e^iWF^fmi +\\G{t)f{tW u )dt 

Jo (2.7) 

<cb,T)ii/r Lr(0oo; ^ 8)) . 

The obtained estimate (|2.7[) pertaining to the integral on (0,T) is now used 
to derive the following one on (T, 00): 



/ e ^||S*e^ t /Wll^di = E / e^ t \\B*e A ' t f(t)f u dt = 

J T k=1 JkT 

00 „t 

= Ve SptT / e 5pT || J B*e AV e A * feT /(r + feT)||P dr 

fc=l ^ 

OC 

<c( P ,T)^e^||e^ feT /(. + fcT)||^ (0co; ^, E)) 

k=l 

00 

< c(p,T)M^^- 5)pkT \\f\\l^oo-MA^) ■ ( 2 - 



k=l 
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Since S < u, the series in (|2.8p is convergent, and there exists a constant C p , 
which depends only on p, such that 



e^We^fmidt < C p \\f\\l ri0tOo . MA , n) . (2.9) 



T 



Assume now on p and r the constraints arisen so far during the proof, namely 

1 p 

1 < p < — , < r < oo . 

7 1 — jp 

The same arguments used to find (I2.6[) provide 

/■oo 

e 5pt ||F(i)/(i)||£ ctt<M p / t-^e-^- S)pt \\f{t)\\ p Y dt 



<c( P )\\f\\ p Lr< 



IT 

I 

which in view of (|2.9I) implies as well 



/•OO 

e^\\G{t)f(t)\\ P udt= / e ^||[F(i)-iTO/(*)ll^ 
Jo 

< c (P. r )ll/llz,-(0,oo;27(X-)) ' 

The above estimate and (|2.4[) finally establish (|2.3[) . 

The proof of (|2.3p in the case r = +oo is similar (even simpler), hence it is 
omitted. □ 

2.2 Regularity of the input-to-state map and its adjoint 

A fundamental prerequisite for all the computations performed in the paper is a 
detailed description of the regularity properties of the 'input-to-state' mapping, 
that is the mapping L which associates to any control function u(-) the solution 
to the state equation (jl.ll) with yo = 0. Namely, L is defined as follows: 



(t) := f e (t - s)A Bv{s)ds 
Jo 



t rt 



F(t- s)*v(s)ds+ I G{t- s)*v{s)ds ^ 2 ' 10 ^ 



L : v i-> Lv , Lv 



o Jo 
= L (1) v(t)+L {2) v(t), t>0. 

We begin by recalling the main regularity properties of the integral operator 
L(iy We note that in view of the key estimate (|1.6p satisfied by the component 
F(t), the proof of the statements of Proposition 12.41 below follows a pretty 
standard route. In fact, it employs the same arguments used in the study of the 
input-to-state map pertaining to parabolic-like dynamics or systems yielding 
singular estimates; see [30] and [28] . 

Proposition 2.4. Let L^j be the operator defined in (|2.10p . Then, the following 
regularity properties hold. 

(i) -L(!) maps continuously L 1 (0,cxd; U) into L r (0, oo;Y) for each r £ [1, — [; 
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(ii) maps continuously L p (0, oo; £7) into L r (0,oo;Y) for each p 

e [p, x _ (1 p _ 7)p ]; 

(iii) maps continuously L 1 ^ (0, oo; [/) into L r (0, oo; y") /or eac/i r E [j^z> 



7 ' 

(iv) Lfys maps continuously L p (0,oo]U) into L r (0, oo;F) n Cb([0, oof; Y) for 
each p E] , oo] and r E [p, oo] . 

Proof. First, an easy application of Holder inequality and Tonelli Theorem 
shows that Lm maps continuously L p (0,oo;U) into L p (0, oo; F) for each p E 
[1, oo]. Next, property (i) follows directly by Holder inequality; property (ii) is a 
consequence of [20j Theorem 383] and interpolation, and finally properties (iii) 
and (iv) follow again by the Holder inequality and interpolation. □ 

The analysis of the operator L(2) is more tricky. It exploits the distinct 
regularity properties of the component G(t) pointed out in Proposition 12.31 

Proposition 2.5. The following properties hold true: 

(i) L(2) maps continuously L 1 (0, oo;[/) into £ s (0, oo; [D(A* e )]') for each s E 

[1,1/7); 

(ii) L(2) maps continuously L p (0, oo;C/) into L r (0,oo;Y) D C{,([0,oo[;Y) for 

each p E (1, oo] and r E [p, oo]. 

Proof, (i) The proof is based on a duality argument. Let u E L 1 (0, oo,C/) 
be given and let w E £ r (0, oo; T>(A* e )), with the summability exponent r to be 
choosen appropriately later. We take the duality pairing ((L(2)u)(t), w(t))-p(A*')i 
integrate on (0, oo) and rewrite as follows: 

POO rt 

(L i2) u{t),w{t)) v{A ,. } dt = ( G(t~T)*u(T)dT,w(t)) v{Ate) dt 

OO pOC 

(u(r),G(t - T)w(t)) (J dtdT 

r-OO 

(u(r), / G(a)w(a + t) da) jj dr , 
jo Jt ' 

which gives 

{L {2)U {t) 1 w{t))- D(A ,e ) dt < ||w||ii(o,oo;cn sup / \\G(a)w(a + t)\\u da 

T>Q l Jt 

(2-11) 

We then proceed with the estimate of the integral on the right hand side of 
(|2.11[) . We insert the exponential weight e St , apply the Holder inequality first 
{p 1 denotes the conjugate exponent of p) and utilize the estimate (|2.3p of Propo- 



14 



sup 

>0 



i/p 



sition l2~3"l next. to find 

(L(2)u(t),w(t)) V ( A **) dt 

< ||«|Ui(o,ooitr) (J e~ Sp ' a dn) 1/P s- 

< Ct ||m||l1(0.oo;(7) SUp ||li?||l,f ( Ti0o; x)M»«)) 

T>0 

= Ct ||«[|ii(0,oo;l7) IMI-L^O.oo:©^")) • 

(2.12) 

Note carefully that Proposition 12.31 applies with any p S [1,1/7) an d to S 
L r (r, 00; X>(A* C )), provided the summability exponent r fulfils the constraint 
Therefore, (I2.12j) implies by duality 



L (2) G £(i x (0, 00; C/), i s (0, co; [P(A* £ )]')) , 
where s is the conjugate exponent of r and the estimate 

\\L(2)u\\l°(0,ooiIV{A-")]')) < C [[tt||ii(0,oo;C/) ) 

holds with a constant C depending only on p and T. Thus, 

P 



readily implies 



1 < s < 



while it is elementary to check that 



p- (l-7p) ' 



g 1 
max = — 

pe[i,^) P~ (1-7P) 7 



which confirms the constraint s G [1, 1/7) and thus completes the proof of part 

(i) . 

(ii) For fixed p £]l,oo[, using Proposition 12. II and Hypothesis ll.4f ii). we com- 
pute for any <p g C£°(]0, 00 [; Y): 



{L {2) v{i),y{tj) Y dt 



Jo 



(v(s),G(t - s)ip(t)) u dsdt 



< 



< 

<c(p) 



|e- tfCt - ) «(-)[[w(o I t i col|e J -G(.)¥'(t)[[ L p'(o,t ! EO 



< 



«(0llLP(O.t:C/)ll t P(i)l|y 



From here, using Holder inequality and Tonelli Theorem, it follows easily that 



as well as 



(L (2) v(t),¥>(*)) y * 



(£ (2) t;(t),v»(*)) y ^ 



< IHIi,p(0,°o;[/)IMIlp'(0,oo;Y)> 



< \\v \\l»(0,oo;U)\\ 1 P\\l~>-(0,oo;Y)> 



which proves, by density, that Li 2 \V € L p (0, 00; f7) n i°°(0, 00; U). In addition, 
the argument of [3J Proposition B.2] yields continuity. The result then follows 
by interpolation. □ 
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The major regularity properties of L are collected in the following Proposi- 



tion. 

Proposition 2.6. The operator L enjoys the following properties: 

(i) L maps continuously -^(O, oo; U) into L r (Q, oo; [D(A* e )]') for each r G [1, — [; 

(ii) L maps continuously LP (0, oo; U) into L r (0,oo;Y) for each p e]1,j^t ( [ and 

r e b) i-(i P - 7 ) p ]/ 

(iii) L maps continuously (Q, oo; U) into L r (0, oo; Y) for each rGfy^— , oof; 

(iv) L maps continuously L p (0, oo; U) into L r (Q, oo; Y) nC(,([0, oof; Y) for each 
p e] y^— , oo] and r 6 [p, oo] . 

(v) L maps continuously L r (0, oo; £/) tuit/i r 6 [g', oo] into C(,([0, oo[; 2?(A e )). 

Proof. Parts (i)-(ii)-(iii)-(iv) readily follow combining the results established in 
Propositions 12.41 and 12.51 

In order to prove (v), let r > q' and let u £ L r (0, oo; U) be given. Then, by 
Proposition 12.21 for each z S £>(A* C ) and i > we have: 



|(Lw(i),A* e z)y| 



(u(t-s),B*e A * s A* c z) c/ ds 



(e~ ds u(t - s),e 5s B*e A * s A* 6 z)u ds 



< 



< c||e- 5 u(t-.)llL 9 '(o,* ;C /)l|e 5 ' B * eA *'^^ll^(o^) < 

< c(g,r)||u|| Lr(0iOO;C/) ||z||y , 

so that Lu G L°°(0, oo; £>(A e )). Now if t > r > we also have 



\{Lu(t) - Lu(T),A* e z)^ 



{u{t~ s),B*e A ' s A* e z)uds- I (u(t - s), B*e A * s A* e z)u ds 



(e~ Ss u(t - s), e Ss B*e A * s A* e z)u ds 



e- 5s (u(t - s) - u(t - s), e fls B*e A s A* E z) c/ 



< 



< C 



(g,r) (lie" 5 u(t - -)\\ Lq > {T>tm + \\e- S \u(t - •) - u(t 



IL?'(0,t;C7) 



||e e -B*e A *-^ e z|U, (Ojtit0 



< 



< c(g,r) (||w|Ur( 0it _ r ;[/) + ||it(- + i-r) - u(-)IU<-(o,oo;i7)) ||z||y, 
and the result readily follows. 



□ 



We recall that the adjoint L* of the operator L is defined as follows: 

/>oo 

L*:v^L*v, L*v(s)= B* e A * ^ v{t) dt, s>0. (2.13) 

J s 

As it is well known, the regularity analysis of L* is also central to the study of 
the optimal control problem. 
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Proposition 2.7. For the operator L* defined in (|2.13[) the following properties 
hold true. 



(i) L* maps continuously L x (0, oo; Y) into L r (0, oo; U) for each r G [1, -[; 

(ii) L* maps continuously L p (0,oo;Y) into L r (0,oo;U) for each p 

andr G [p, x _ (1 p _ 7)p ]; 

(iii) L* maps continuously L p (0, oo; y) into £ r (0, oo; {/) /or eac/i p € [-j^-, oo[ 
and r G [p, oo[. 

(iv) L* maps continuously L p (0, oo; T>(A* e )) into L°°(0, oo;U) for any p G 
] 137.0°]; 

(v) L* maps continuously L 1 (0, 00; [D(A C )]') into L r (0,oo;?7) /or a/Z r G [l,o]- 
Proof. It follows by Proposition 12 . 61 by duality: indeed, we have for each p > 1 



{L*v(t),u(t))u dt 



(v(t),Lu(t))u dt 



Lp(0,oo-U)\\ Lu \\lp'(0,oo-U) 



Now to prove (i) we take p = 1, so that p' = 00, and by Proposition I2.6f iv) we 
get 

f°° / 1 i 

f {L*v{t),u{t))udt < c||u[[ i i( 0)00 ;i7)[|ii|| i r'(o i01o; t/) Vr'el- .°° . 

T 

which means L*v G £ r (0, 00; {/) for each r G [1, 1/7). 

Similarly, to prove (ii) we take p G (1, 737), so that p' G (—,00): then, since 
P' = i-(i- 7 )r' if and only if r ' = i+(i-7)p' 1 Pr °P° sition IHIii) y ields 



(L*v(t),u(t))udt < ||u||iP(0,cc;t/)l|w||ir'(0 1 oo;I7) Vr ' G 



1 + (1 - 7 )p' 



7'f 



which means L*w S £ r (0, 00; J7) for each r £ [p, 1 _( 1 °_ 7 - )p ]- 

To prove (iii), we take p G [3^37,00), so that p' G (1, -]; since again p' 

i_(i"_ 7 ) r ' ^ an d om y if r ' = i + (i'_ 7 )p' i by Proposition I2.6f i) we obtain 



(L*v(t),u(t))udt < c||w|| L p (0 , oo;C/) ||m|| l ^ (0 oo . [/) Vr'G 



l + (l-7)p 



7'P 



which means L*w G £ r (0, 00; {/) for each r G [p, 1 _ ( - 1 °_ 7 - ) p ]- 

The assertion (iv) is the dual statement of (i) of Proposition 12.61 

Finally, (v) follows again from the assertion (v) of Proposition l2~6l bv duality. 
We take v G L 1 (0,oo;[V(A t )]'), which by Remark A.l(ii)] means v{t) = 
A* e w(t), with w G L 1 (0, 00; y). We then compute 



(L*v(t),u(t))udt= / (A* e w(t),Lu(t)) Y dt 



(w(t),A £ Lu(t)) Y dt < 11^11x1(0,00;^) ||iu||c„([0,oo[;I5(A«)) 
< Cr 1 1 A* ~*V II Li(0,oo;y) |K|U'-(0,oo;i7) > r G W, °°] . 
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where in the latter estimate we used (v) of Proposition 12.61 The above shows 
that L* maps i 1 (0, oo; [P(A e )]') continuously in the dual space of L r (0, oo; U), 
for any r G [g',oo], that is U (0, oo; U) with r' G [l,q], thus concluding the 
proof. □ 

Remark 2.8. All the regularity results provided by the statements contained in 
Proposition 12.61 and Proposition 12.71 may be easily extended to natural analogs 
in L^-spaces, if <5 G [0,w A 77) is given, maintaining the respective summability 
exponents p. The proof is omitted. 

In the following Proposition we collect a couple of regularity results which — 
in view of Remark 12.81 — are in essence contained in assertion (v) of Proposi- 
tion 12.61 However, since (|2.14p below will play a crucial role in the proof of 
well-posedness of the Algebraic Riccati equations, its statement is given explic- 
itly, along with a distinct proof. 

Proposition 2.9. For any 8 G [0,wA r/[, the following regularity results are 
valid: 



R*RL G C(L q s (0, 00; U),Lf(Q, 00; V{A* e ))) ; 
L*R*R G £(Lj(0, 00; [D(A* e )]'), L|(0, 00; U)) . 

Proof. Let u G L q s ' (0, 00; U) and z G T>(A e ). We have 

\(R*RLu{t),A e z) Y \ = f {u{a),B*e A " it - a) R*RA e z) u da 
Jo 



(2.14a) 
(2.14b) 



< 



< Ce 



.(tr), e Sit ^ ) B*e A " (t - a '>R*RA e 



< 



dt < 



-St ||„<5- 



e m IIl<!'(o,oo;;7) 



e 5 B*e A *R*RA e z\ 



Ll(0,oo;U) 



< 



-St 



<Ce ~ ||e" u|| i9 ' (0)Oo;£/) 
where 5 G]0, a; A n[ and in the last estimate we used Proposition ^. 21 Therefore, 

\(e 5t R*RLu(t), A*z) Y \ < C \\e s u\\ Lq , {0 oo;U) \\z\\ Y , 

i.e. A* e e 5t R*RLu G L°°(0, 00; U) and 

\\R*RLu(-)\\ Lr (a^viA**)) < c\\u\\ L t> (0)OO;tr) • (2-15) 
The statement in (|2.14b|l follows from (I2.14ap by duality. 

□ 



3 The optimal state semigroup, the optimal cost 
operator 

3.1 The optimal pair 

Since we are dealing with a classical linear-quadratic problem, the existence of 
a unique optimal control minimizing the cost functional (|1.3|) follows by convex 
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optimization arguments. In addition, the Lagrange multipliers method yields 
the optimality condition for the optimal pair (y, u) G L 2 (0, oo; Y) x L 2 (Q, oo; U), 
that is 

u = -L*R*Ry, (3.1) 

where 

y = e A -y + Lu. (3.2) 

Then, owing to the regularity provided by Propositions 12.61 and 12.71 as well as 
to the decay assumption on the semigroup e At , we can appeal to a classical 
bootstrap process: we set 

po = 2, p n+1 = y~ — - — , < n < N , 

1 - (1 - j)p n 

where N G N is the first positive integer such that pn > j^—', such an integer 
does exist because 

(1- 7 K . 4(l- 7 ) 

Pn+l ~Pn = Pn Z JZ : > — — > . 

1 - (1 - "f)p n 27-1 

Since u G L Po (0,oo; U), by (pT2"j) we know y G L Pl (0, oo; Y), so that we deduce 
as a first step u(-,s;x) G L r (0,oo;[/) for all r < pi as well as, by (|3.1[) . y G 
L s (0, oo; Y) for all s < p%. Thus, after n steps, we get 

u G L r (0, oo; U) Vr<p 2n , y G L s (0, oo; Y) Vs < p 2n +i ■ 

This procedure stops as soon as 2n or 2n + 1 equals N: indeed, if N = 2n 
we get at the n-th step u G L r (s 7 T : U) for all r < px and y G Cb([0, oo[; F), 
so that in the next step we obtain u G i p (0, oo; U) for all p < oo; if N = 2n + 1 
we find directly, at the (n + l)-th step, u G L p (0, oo;C/) for all p < oo and 
y G Cb([0, oo[; 1"). At any step we also find the corresponding bound in terms 
of ||x||y . Thus, we have proved the following Proposition. 

Proposition 3.1 (Statement SI. of Theorem 1 1.5p . For any x G Y there exists 
a unique optimal control u(-) for problem (jl.l[) - (jl.3[) . The optimal pair (y,u), 
which a priori belongs to L 2 (0, oo; Y) x £ 2 (0, oo; U), further satisfies 

yG J L p (0,oo;F)nCf ) ([0,oo);y), u e L p (Q,oo;Y) VpG[2,oo), 

continuously with respect to x G Y : 

||2/|U*>(0,oo;C/) < c p \\x\\ Y VpG[2,oo); ||j/||c b ([o,oo[ ; y) < c\\x\\ Y , (3.3a) 

\\ u \\lp(o,oo;U) < c p\\ x \W VpG[2,oo). (3.3b) 

Let us denote by $(-)cc the optimal state corresponding to an initial state 
x G Y and to the optimal control u(-): 

$(t)x := y(t) = y(t,x;u) . 

The significant basic properties of the family of operators {$(i)}t>o are briefly 
recorded below. 

1. The estimate on the right in (|3.3ap of Proposition I3TT1 yields ||3?(f)||,c(y) < oo 
uniformly in t > 0. 
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2. It is not difficult the show that y(t,x) possesses the transition property 

y(t + a;x) = y(a;y(t,x)) £ C 6 ([0,oo);Y) xeY,i,CT>0, (3.4) 

which renders &(t) a strongly continuous semigroup on Y; indeed, $(0) = I is 
readily checked, as y(0, a;) = x for any x £ Y. We just note that the proof of 
(|3.4j) follows a standard route: it is based on the identity 

y(t; x) + [LL*R*Ry(-, x)](t) = e At x , (3.5) 

which follows from Lu = —LL*R*Ry (this, in turn, is a direct consequence of the 
optimality condition (|3.1[l ). combined with p.2[) . (See, e.g., [501 Theorem 6.25.1, 
p. 626] for an outline of the necessary steps, though in the context of min-max 
game theory.) 

Thus, (|3.4D confirms that $(£) is a Co-semigroup on Y, which introduces its 
infinitesimal generator, that is the linear, closed operator Ap : V(Ap) C Y — > Y 
defined by 

<f>(t)x — x 

V(A P ) := \ x £ Y : lim — ^ exists in Y\ , 

^"° + t (3-6) 
,. $(t)x -x ^, A . y J 

A P x := lim — ^ , x £ V(A P ) . 

t^o+ t 

The domain T>{Ap) is by its nature dense. 

3. Since $(-)a; £ L 2 (Q, oo;Y) and is a Co-semigroup, then according to 
Datko's Theorem = e Apt is further exponentially stable, which confirms 

CUD. 

In summary, we have the following 

Proposition 3.2 (Statement S2. of Theorem 11.51) . The optimal state y(t,x) 
defines a strongly continuos semigroup e Apt in Y , t > 0, which is exponentially 
stable. 

A deeper insight into the additional differential properties of the optimal 
state semigroup <!>(£) = e Apt — the ones which will ultimately allow us to estab- 
lish well-posedness of the ARE — is given in Section @] 

3.2 The Riccati operator P 

The Riccati operator P is initially introduced, as usual, explicitly in terms of 
the optimal state semigroup <&(f): 

y* oo 

Px:= e A ' t R*R<S>{t)xdt, x£Y. (3.7) 
Jo 

We record first the basic properties of the operator P. 

Proposition 3.3. The operator P belongs to C{Y) n C(V(A £ ), V(A* e )). 

Proof. With x £ Y we just utilize (11.21) of Hypothesis 11.11 and the continuity 
property (|3.3ap of Proposition 13.11 to find 

/>oo 

ll^llr < /_ \\e Att \\c(Y)\\R*R\\c(Y)\\Ht)x\\ Y dt< 
i Mce~ ut dt \\x\\ Y <c'||x||y. 
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Next, for x G V(A e ) we write $(i) = e At x+Lu(t) and employ the assumption 
(ll.8[) on the observation operator i?, the regularity provided by (|2.15[) and (I3.3a[) 
to obtain 

/•OO 

\\Px\\ v{A , n < / \\e AH \\ c{Y) \\R*R[e M x + Lu{t)]\\ v(A , t) dt 
Jo 

/>oo 

/ Me- ut \\A* e R*RA- £ \\ r(Y ,\\e At A e x\\ Y dt 



\ C (Y)\\e ax\\ y < 

+ / Me- w *||irilLu||Lco ( o l00i i,( J i«))dt 
Jo 

< C [||x||z>(A«) + INIz,«'(0,oo;£/)] ^ c \\ x \\v(A') ■ 



□ 



By the very definition of P in (|3.7[) it easily follows — using the optimality 
condition (13.11) — that the optimal cost is a quadratic form on Y, that is (-Pa;, x); 
this motivates the term optimal cost operator used for P. The proof of this fact 
is briefly recorded below. 

We simply compute, for any x G Y, the cost corresponding to the optimal 
pair, that is 

J{u)= j ^{Ry{t,x),Ry{t,x)) Y +(u{t,x),u{t,x)) u ^dt 

p OO 

(Ry(t,x),Ry{t,x)) dt- / (L*R*Ry(- ,x)(t),u(t,x)) dt 

Jo 

/•OO 

(R$(t)x, R[e At x + (Lu)(t)]) dt - (R<f>(t)x, R(Lu)(t)) Y dt 

Jo 

= (R<S>(t)x,Re At x) Y dt= (J e A ' t R*R$(t)x dt, x)^ . 

Owing to (|3.7p . the last equality shows that 

J(u) = (Px,x) Y VxeF. (3.8) 

Remark 3.4. In much the same way one gets J(u) = (x, Px)y, which combined 
with (|3.8p shows that P is a self-adjoint operator; moreover, since J{u) > for 
any x G Y, (|3.8|) establishes as well that P is non-negative. 

Similarly, one may show that 
(Pxi,x 2 ) Y = J {(Ry(t,xi),Ry(t,x 2 )) Y + (u{t,xi),u(t,x 2 )) u } dt, xi,x 2 eY. 

We now introduce the function p(t, x) defined by 

/OO 
e A * (T - t) R*Ry(T,x)dT t > , x G Y , (3.9) 

which is the unique solution to the Cauchy problem 

(p'(t;x) = -A*p(t)-R*Ry(t,x) 

\ lim p(t ;a: )=0. (3 - 10) 
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A classical and elementary computation shows that if we rewrite P in (13. 7|) 
as a time-dependent function, that is 

/oo 

and next apply the above formula when x is replaced by <&(t)x, namely: 

/oo y* oo 

e^*(T-*) j R*ij$( T _ t )$( t ) a; dt = J e A * {T ~ t} R*R<P{T)xdt, 

then we immediately see that 

P$(t)x=p(t;x) Vt>0,x£Y, (3.11) 

which establishes a relation between the Riccati operator P and the function p. 

The above formula p. lip is the starting point in order to derive the feedback 
representation of the optimal control — a property which is central to solvability 
of Problem 1 1.3 1 Indeed, if we return to the optimality condition (|3.ip . and write 
explicitly the integral operator L* , we know that 

/oo 
B*e A " {T - t) R*R<I>(T)xdT Vx E Y . (3.12) 

Thus, a formal computation yields 

/oo 
e A*{r-t) R * m ^ xdT= _B*p(t;x), (3.13) 

which combined with p. lip would imply 

u(t,x) = —B* P$(t)x for any x G Y and for a.e. t > 0. 

However, going from (|3.12p to (|3.13p necessitates a deeper technical justification, 
which is found in the next section. 

3.3 The gain operator B*P 

The technical issue raised in the previous section — namely, the need for a rig- 
orous justification of the first equality in (|3.13p — as well as our next and major 
task (that is to show that the optimal cost operator P does satisfy the algebraic 
Riccati equation corresponding to the optimal control Problcm ll.3p . require that 
we are able to give a proper meaning to the gain operator B*P. We accomplish 
this goal by introducing a linear operator, which will eventually coincide with 
B*P, that is shown to be bounded on a dense subset of Y (and yet unbounded 
on Y). The regularity result for the operator R*RL set forth in Proposition ^. 91 
along with Proposition 12 . 31 provide the tools. 



Theorem 3.5 (Statement S4. of Theorem 1 1 . 5|> - Let e be such that Hypothe- 
sis \1.4^ iii) holds true. Then, the following statements are valid. 

(i) The integral 

poo 

Tx:= B*e A ' t R*R$(t)xdt (3.14) 
Jo 

defines a (linear) bounded operator from T>(A e ) into U. 
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(ii) The operator B*P 7 formally defined by 

B*Px:=B* e A ' t R*R${t)xdt 1 
Jo 

coincides with T on T>{A e ), and hence 

B*P e C(V(A e ),U). (3.15) 

Proof, (i) Let x € V(A e ). According to the decomposition (|1.5[) for B*e t , the 
integrand in (|3.14p is split as follows: 

B*e AH R*m{t)x = F{t)R*m(i)x + G{t)R*Re At x + G(t)R* RLu(t) , (3.16) 

with a further splitting due to (|3.2[) . We will show that each summand on the 
right hand side of (|3.16p is a function in L 1 (0, 00; U). 

First, the basic Hypothesis II. 4f i) combined with the regularity Q3.3a[) of the 
optimal state readily imply 

\\F(t)R*R$(t)x\\u < Nt~^ e-^\\R*R\\ ciY) \m)x\\ Y < ct^e"* ||x||y , 

which shows 

nOC 

/ \\F(t)R*R^(t)x\\udt <c\\x\\ Y VieF. (3.17) 
Jo 

Next, we recall Proposition O and utilize (123)) with p = 1, f(t) = R*Re At x 
and r = 00, where R*R is subject to Hypothesis 1 1 ,4f hi) (b) . thus obtaining 

poo 

/ \\G(t)R*Re M x\\udt < c sup \\R*Re At x\\ v{A „) < 
Jo t>0 

< C\\A**R*RA*-*\\ c{Y) \\x\\ v{A < ) , 

so that 

/>OC 

/ \\G{t)R*Re At x\\ u dt<c\\x\\ T){At) Vx e V{A e ) . (3.18) 
Jo 

Finally, for the third summand in (|3.16p we apply once more Proposition ^. 31 
this time with f{t) = R*RLu{t). Notice that the membership R*RLu(-) G 
L°°(0, 00; T>(A* € )) follows from (|2.14a|) . in view of the regularity established in 
Proposition 13. II which provide as well the appropriate estimate. Therefore, 

pOO 

/ \\G(t)R*RLu(t)\\u dt<c sup \\R*RLu(t)\\ v{A , t) <c\\u\\ Lq , (0 u): 
Jo t>0 

which yields, in view of (|3.3b[) . 

poo 

/ \\G(t)R*RLu(t)\\udt<c\\x\\ Y VieY. (3.19) 
Jo 

The obtained estimates ([5Tf| . (|3~T51) and (l3~lT))) show that 



\Tx\\u = 



/ B*e A ' t R*R<S>(t)xdt < c ||a?||x» { A-) Va; E T>(A e ) . 
Jo u 
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that is T G C(V{A £ ),U). 

(ii) We return to the decomposition (11.51) , where by Hypotheses ll.4f i) and (iii) (a) 
we know 

F(t)€C(Y,U), G(t) G C(V(A* e ),U) , 
respectively. On the other hand, we also have for t > 

B*e A ' t x=[e At B]*x Vx G V{A*) , (3.20) 

where 

[e At B]* eC(V(A* e ),U); 

see [H Lemma A.2]. Thus, the operators F(t)+G(t) and [e At B]* are well defined 
on V(A* e ) and both coincide with B*e A " 1 on 2?(A*). Hence, if x G X>(v4 e ) we 
can write for all v G 17 



{Tx,v)u 



= / ([e jit S]* J R*i?$(t)ar,«)i 7 dt = 
Jo 

/•oo 

= / (i?* J R$(i)x,e Ai J Bw) I , ( ^»» ) ,p (A , e)] ,^ = 
Jo 

/•OO 

= lim / (e A ' s R*R^t)x ) e M Bv) v{A ** uv{A ,c ) ydt 



-lim / (e A t R*R^{t)x,e As Bv) v{A , %[D(A , c) ydt = 
s ^ 0+ Jo 

= lim (Px,e As Bv) v(A , t)AV{A , t)V = lim ([e As B]*Px, v) u (3.21) 

s— >0+ s— »0+ 

=: (B*Pa;,v)y. 

This shows that for all x G V(A e ) the operator _B*P is uniquely defined as a 
weak limit in [/, and coincides with 7", i.e. B*Px — Tx. By part (i), 7?*P = 
T E C(V(A e ), U) which concludes the proof of (ii). □ 

Remark 3.6. To give a deeper insight into the previous result, we summarize 
the major steps of its proof, complemented by few additional technical remarks. 
First, the operator T defined by (|3.14[) is shown to be a bounded operator from 
V{A e ) in U. Next, the equivalence between (Tx,v)u and the limit in (|3.21[) 
(for any x G T>(A e ) and any v G U) shows that Tx coincides with the weak 
limit w-\im s ^, Q +[e As B]* Px in U. Notice that the operator [e As i3]* is bounded 
from D(A* e ) in U, and yet its norm blows up, as s — > + , like s -7 ; see 
Lemma A.2]. 

Finally, the motivation for denoting by B*Px the aforesaid weak limit is the 
following. When Px belongs to V(A*) rather than just to T>(A* e ), we know 
that [e As B]*Px = B*e A " s Px -> B*Px, because e A " s Px -> Px in V(A*) as 
s — > + , whereas by its very definition B* G L{V(A*),U). We note that the 
membership Px G D{A*) does hold, at least for x G T>(Ap), as we will see later 
in Lemma 14.151 

3.4 The feedback representation of the optimal control 

On the basis of the analysis carried out in the previous section, we would intend 
now to derive the feedback representation of the optimal control. The obtained 
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property (|3.15[) for the gain operator suggests that we preliminarly assume x G 
T>(A e ) and show the validity of the feedback formula for these smoother initial 
data. In fact, with x 6 V(A e ), we can exploit the additional regularity of the 
optimal state, which is made clear in the following Proposition, also of intrinsic 
interest; the proof is postponed to Appendix |A"1 

Proposition 3.7 (Statement S8. of Theorem II. 5|) . If x G V(A t ), the optimal 
state &(t)x belongs to T>{A IL ) for all t > 0. More precisely, 

x G V{A e ) $(• )x G L p {0, oo; V{A e )) n C 6 ([0, oo); T>{A e )) (3.22) 

for all p G [l,oo], continuously with respect to x. Consequently, 

x G D(A e ) ==> u G Cf,([0, oo); U) . (3.23) 

Assuming x G V(A e ), by Proposition[3T7]we know that $(• )x G £ p (0, oo; £>(A £ )) 
and we can appeal to the same arguments used in the proof of Theorem 13.51 to 
give a rigorous justification to p. 131) . It is now true that (I3.11[) . combined with 
(|3.13p . provides the feedback representation of the optimal control, initially for 
x G V(A e ): 

u(t,x) = -B*P$(t)x for any x G V(A e ) and for a.e. t > 0. 

This formula is easily extended to all of x G Y, as shown in the Proposition 
below. 

Proposition 3.8 (Statement S9. of Theorem II .5[) . The (pointwise in time) 
feedback representation of the optimal control 

u(t, x) = -B*P$(t)x , for a.e. t> 0, (3.24) 

is valid for any initial state x G Y. 

Proof. We know that the relation (|3.24j) holds for a.e. t > and for any x G 
T?(A e ). Thus, recalling the (continuous with respect to x) estimate (|3.3b|) for 
the optimal control, we obtain 

\\B*P^(-)\\ LP(0<oo . u) = ||£(-)|U*(o,oo;tO < C\\x\\ Y VxGP(A e ). (3.25) 

By density, the operator B*PQ( ) is extended to a bounded operator from Y to 
L p (0, oo;U). Consequently, 

u(t,x) = —B* P&(t)x for any x G Y and for a.e. t > 0. 

i.e. the feedback formula (I3.24|) is extended as well to all x G Y. □ 

4 Towards well-posedness of the Algebraic Ric- 
cati Equation 

In this section we show that the optimal cost operator P introduced in Section[3] 
does solve the ARE (14.45[) corresponding to the optimal control Problem 11.31 
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Owing to formula (|3.7p , this issue is strongly related to certain differential prop- 
erties of the optimal state semigroup $(i), which are discussed in Section I4T21 
culminating with the statement of Corollary 14. 141 

We will see that a major challenge arises on the operator-theoretic side, as the 
mapping A = I + L*R*RL — which, like in previous theories, is an isomorphism 
on L 2 (0, oo;U) — fails to be an isomorphism in L q (0, oo;U) or in L'(0,oo;£/), 
whilst it is required to admit a bounded inverse acting at least on the latter 
space. A solution to this question is given by the distinct result established 
with Theorem 14. 121 based on Lemma T4. 101 

To begin with, we provide a preliminary description of the the generator Ap 
of the optimal state semigroup $(t). 

4.1 The optimal state generator A P . Basic facts 

Let Ap : T>{Ap) c Y — > Y be the (optimal state) generator defined by (|3.6p . 
i.e. the infinitesimal generator of the strongly continuous semigroup $(£). With 
x G T>(Ap), we know that 

3 Urn ( *® X ~ X ,z) =: (A P x,z) Vz e Y , 

which — by using the representation of the optimal state &(t)x in terms of the 
optimal control u(t, x) — is readily rewritten as follows 

3 lim ( eMx ~ x + 1 f e A l t - r )Bu(T-,x)dT t z)=(A P x,z) VzeY. 

i^0+ V t t J J 

Still with x G T>(Ap) and taking now z G T>(A*), we find 
-($(t)x - x,z) Y = 

= j([e At - QA^x, A*z) Y + \{J e^-^A^Buir, x) dr, A*z) y , 
that is, by setting w = A*z, 

~mt)x-x,(A*)- 1 w) Y = 

1 1/ r* ^ 

= _([ e ^* _ I]A- l x, w)y + j[J A' 1 Bu(t, x) dr, wj ^ . 

Thus, if we let t — » + in (|4.1D . we see that the left hand side tends to 
{Apx, (A*)~ 1 w)y , whilst the first summand in the right hand side converges 
to (x,w)y- As for the second summand in (|4.1j) . we employ the feedback rep- 
resentation of the optimal control (|3.24l) and observe that x G T>(Ap) yields 
<Z>(-)x G C([0,T],V(A P )), which implies P$(> G C([0, T],V(A*)), since P G 
£(T>(Ap),V(A*)) (this fact will be shown later: see statement (i) of Lemma l4.15[) . 
Consequently, B*P$(-)x G C([Q,T],U) and the map 

r i-> e Al 'A- x Bu{t, x) = -e Ar A~ X B B*P<S>{t)x (4.2) 
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is continuous on [0, oo) with values in Y, thereby ensuring that 

3 lim - / e A{t - T) A~ 1 Bii(T,x)dT = -A- 1 BB*Px, xEV(A P ). 
t^o+ t J 

Now in view of the remarks above, (|4.ip yields for any x E T>(Ap) 
{A P x, (A*)- 1 w) Y = (x, w) Y - {A^B B*Px, w) Y Vw 6 Y , 

that is 

(A^Apx, w) Y = (x - A~ x B B*Px, w) y Vw e Y , 
meaning that x - A~ 1 B B*Px E V(A), with 

v4 P x = A(x - A~ l B B*Px) Vx E P(A P ) . (4.3) 

In addition, semigroup theory provides the basic differential property 

— $(t)x = All - A- 1 BB*P]<f>(t)x = <$>(t)A[I - A- 1 BB*P]x Wx e V(A P ) . 
at 

We have established Statement S5. of Theorem ll.5l which is recorded in the 
following Proposition. 

Proposition 4.1 (Statement S5. of Theorem 1 1.51) . The infinitesimal generator 
Ap of the (optimal state) semigroup $(£) defined in (|1.10[) coincides with the 
operator A(I ~ A~ 1 BB*P), on the domain 

V{A P ) C {x E Y : x - A- 1 BB*Px E V(A)} (4.4) 

i rt 

C {x € Y : 3w-lim t ^ 0+ - / e Mt ~ T) A- 1 BB* P<I>(t)x dr (4.5) 
t Jo 

Remark 4.2. We observe that for x 6 T>(Ap) the weak limit 

lim if f e A ^ T) A- 1 Bu(T,x)dT,w) Vw E Y 
t^o+ t\J Jy 

defines a linear operator, which we denote by T, that is 

1 /"* 

Fx := w-lim t _, + - / e A{t ~ T) A' 1 Bu(t, x) dr (4.6) 
t Jo 

1 f* 

= - w-lim t _ >0+ - / e A(t-T) A -1 .B£*P*(r)a:dT, (4.7) 
t Jo 

with dense domain in Y, as £>(Ap) C T>(T). Then, in the discussion leading 
to the statement of Proposition 14. II we have shown that when x E T>(Ap) then 
Tx = -A~ 1 B B*Px, i.e. T coincides with the operator -A- 1 BB*P on V(A P ). 

We note, in addition, that if x E T>(A e ), owing to (I3.23[) of Proposition 13.71 
the map in (|4.2[) is continuous as well, so that the weak limit (|4.6[) is strong, 
V(A f -) C V(T), and we find again 

Tx = -A- 1 BB*Px x E V(A e ) . 
27 



Remark 4.3. Pretty much in the same way we obtain for x G 'D(A) the exis- 
tence of the weak limit 

lim -( [ A- 1 e A(t - T) Bu(T 1 x)dT,'w) VwGY, 
t^o+ t Wo Jy 

which defines a linear operator, which we denote by Tp, that is 

1 r* 

T P x := w-lim t ^ + - / Ap 1 e A ^ t - T ' Bu(t, x) dr 
t Jo 
i rt 

= - w-lim t ^ 0+ - / A p 1 e A{t - T) BB*P$(T)xdT , 
t Jo 

with dense domain in Y, as T>(A) C V{Tp). Moreover, we have 

{Ax, A* p ~ 1 w)y = (x, w)y — (Tpx, w)y Vio G Y , 

i.e. x - Tpx G V(A p ), with 

Ax = A P (x - T P x) Vx G V(A) . 

It is interesting to note that although an explicit representation of Tpx for 
x G T>(A) is missing, we will later see that the key issue is not so much to be 
able to deal with Tpx by itself, but rather to find the "right" representation 
of the operator e Apt ApTp, along with an appropriate regularity (in time and 
space). Specifically, we will establish 

e Apt A P Tp = e Apt B B* P : V{A) — ► [V{A*J for a.e. t > 0, 

a property which is central to the proof of Corollary 14.141 and then to well- 
posedness of the algebraic Riccati equations. 

Remark 4.4. It is important to emphasize that the statement of Proposition l4~T1 
cannot be improved to assert the equivalence between T>(Ap) and the subset of 
Y in (|4.5[) (which is characterized by the existence of the discussed weak limit 
Tx). In fact, assuming x G Y is such that Tx does exist, we only obtain that 

3 l im ( *® X - X )Z ) y z eV(A*), 

t^o+ V t / 

which is not sufficient to conclude that x G T>(Ap). 

One may also wonder whether the inclusion in (|4.4p becomes an equality, 
possibly adding the restriction x G T>(A e ). However, this is not the case. In 
fact, assuming x G T>(A e ), with x—A~ 1 BB*Px G T>(A), from the decomposition 

( *® X ~ X ,z) = (±[e M -l](x-A- 1 BB*Px),z) + 

+ (l J eA{t ~ s)BB * P \ x - ^(s)x],zj , t > , z G Y , 

we immediately see that the first summand on the right hand side converges to 
A(x — A^ 1 BB* Px) as t — > 0, while to make sure that the second summand is 
infinitesimal we need z G V(A*). Thus, the same conclusion as above follows. 
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Although a full characterization of D(Ap) is missing, the following Propo- 
sition clarifies the relation between T>{Ap) and T>(A e ), as well as between the 
domains of the corresponding adjoint operators. This result will be also em- 
ployed in the proof of Corollary 14.141 

Proposition 4.5. The following inclusions are valid, provided e < 1 — 7; 

T>{A P ) C V(A e ) , V{A* P ) C V{A* e ) . (4.8) 



Proof. We prove first the latter inclusion in (|4.8p , which is in addition of central 
importance in the proof of Corollarv l4.141 the former can be shown using similar 
arguments. The proof is based on a well known characterization of the domains 
of fractional powers A* e in terms of the interpolation spaces (Y, T>(A*)) t: 2- The 
idea is to relate first V(Ap) to one of the interpolation spaces (Y,T>(A*)) a , oc . 
These can be also described as follows: 

(Y,V(A*)) a:00 = {xeY : sup t- a \\e A ' l x - x\\ Y < 00}; (4.9) 

te(o,i] 

see [351 Theorem 1.13.2]. 

We aim to show that V(A* P ) C (Y, V(A*)) a ^ OQ , for a G (e, 1). Lets G V{A* P ). 
Then 

[$(*)* - I]x = 0{t) , t 0+ , 

which implies 

(Mty-i]x, z ) Y = o(t), t^o+, 

for all z E Y. With x E 2?(Ap) fixed and any z G Y, we rewrite 

([$(t)* -I]X,Z) Y = (X, [$(<)- I]*)y 

= {x,[e At - I]z)y + (x, [ e A{t - s) Bu(s,z)ds) Y (4.10) 



T 2 (t) 

and focus on the second summand T2, which in turn splits as follows: 

T 2 (i) = (x,L ( i)u(t))y + (a;,£(2)u(t))y . 

(The operators Lu\ — resulting from the splitting of B*e A 1 — have been intro- 
duced in (|2.10p : above, we set L^u{t) in place of {Lu\u{-, z))(t), i = 1,2, just 
for conciseness.) It is readily seen that 

\{x,L {1) u(t)) Y \<\\x\\Y jj—^-\\u(s,z)\\uds 

ft 1 1 \Vp' 
< Ci ||x||y ||w|| L p (0:OO;;7) y ^ _ ^ ip , ds j 

<C x ^-^\\x\\ Y \\z\\ Y . (4.11) 

To achieve the above estimate we have used Assumption 1 1 ,4f i) . the Holder in- 
equality, as well as the continuity property (|3.3b|) established in Proposition ^. II 
Notice that by Proposition 13.11 v can be taken arbitrarily large: in particular, 
here p > 1/(1 — 7) is required, in order to ensure that the exponent l/p' — 7 
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is positive. In addition, as by assumption 1 — 7 — e > 0, then we may choose 
P > (1 — 7 — so that 1/p' — 7 > e, and (|4.12l) reads as 

\(x,L (1) u(t)) Y \ < d i Ql ||z|| y (4.12) 

with Oi\ > £■ 

As for the summand (x, L( 2 )u(i))y , by Proposition I2.5f ii) we know that 
t n- L(2)u{t) is a continuous function on the whole half-line [0,oo), as u E 
L p (0, oo;?7) for any finite p > 1. On the other hand, since we aim here to 
obtain an asympotic estimate of L(o\u(t) as t — > + , we may set t < T. Taking 
the inner product with any x E Y, we find 

\(L( 2 )u(t),x) Y \ = / (G(t- s)*u(s),x) y ds < / |(u(s),G(t - s)x)u\ ds 
Jo Jo 

(4.13) 

< \\u\\lHo,T;U) (J \\G(t-s)x\\ q u d s y /q (jf ld S ) lA (4.14) 

< ||fi||i»(o,ooju) ||x||y l|G(0lk(y,^(o ) T;t/)) i 1/r < ||z||y||x|| y Vx E F , 

where to go from (|4.13[) to (14.141) we applied the Holder inequality with l/p + 
1/q+l/r = 1, and by Assumption ll.4f ii) the summability exponent q, like 
p, can be chosen freely as well. Notice that this makes it possible to render 
1/r = 1 — 1/p — l/q arbitrarily close to 1. The above computations yield the 
pointwise estimate 

\(x,L (2) u(t)) Y \ <C 2 t a *\\x\\ Y \\z\\ Y (4.15) 

with arbitrary a,% < 1. Thus, combining (|4.12l) with (|4.15[) we find that there 
exists a constant C such that 

\T»(t)\ < \(x,L (1) u(t))Y\ + \{x,L {2) u{t)) Y \ 

< C t™^^\\x\\ Y \\z\\ Y = 0(t a )\\x\\ Y \\z\\ Y , 0<t<T, 

with e < a < 1. 

Returning to (|4.10[) . we have so far shown that 

\\{e AH - I)x\\ Y = 0(t) - 0(t a ) = 0(t a ) , i^0+, (4.17) 

which in view of (|4.9[) establishes the membership x E (Y, T>(A*)) at00 for all 
a E (e,l). 

Thus, if we recall from [351 Theorem 1.3.3] the inclusions 

(X,Y) aA c {X,Y) atP c {X,Y) atOB E {X,Y) e ,i , 

which hold for all a, 9, p such that < a < 9 < 1 and 1 < p < 00, we 
immediately conclude that there exists 9 E (e, 1) such that 

xe(Y,V(A*))e,2 = V{A* e )- 

see, e.g., [30l § 0]. Consequently, x E V(A* e ) which shows 2?(A£) C V(A* t ), 
thus concluding the proof. □ 
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4.2 A distinct regularity result pertaining to e Apt 

We discuss here a distinct regularity property of the optimal state semigroup 
<!>(£) = e tAp which will play a crucial role in the proof of well-posedness of 
the Algebraic Riccati Equations (ARE) corresponding to the optimal control 
problem. It is indeed the regularity result established in Proposition 14. 1 31 below 
which will make it possible to differentiate strongly the semigroup e Apt on T>(A), 
as needed to obtain that the optimal cost operator P does satisfy the ARE on 
V(A). 

Remark 4.6. The same issue, that is strong differentiability of e Apt on T>(A), 
was addressed as well in the study of the infinite horizon LQ-problem for abstract 
control systems which yield singular estimates] see |28| Section 3.3] . In that work 
the sought property was easily established in view of the following crucial fact: 
$(t)_B inherited the same singular estimate as e At B. 

To begin with, let us preliminary state the intrinsic regularity of the map 
t M> e At B for the class of control systems under investigation. 

Lemma 4.7. Consider, for t > 0, the operator e At B, defined in U and taking 
values — a priori — in [D(A*)]' . For any S G [0,wAi)[ we have 

e s e A B G £{U : L s (0,oo;[D(A* e )}')) Vs e [1, -) . (4.18) 

7 

Proof. We use a duality argument. With u € U and / E L r (0, oo; V(A* t )), 
1/(1 — 7) < r < 00, we estimate 

) /.oo 

(e st e M Bu,f(t)) v{A „ ) dt = / (e' et u,e^ t B*e A ' t f(t)) u dt 

Jo 

/ f°° \ i/p' 

< (/ e-^dt) \\u\\ u || eV+ e >B*e A -f(-) \\ L * {0 ,oo;U) (4.19) 

< C||u|| ff ||/|| Z r (0lOO; 2, ( A«)), (4.20) 

where 9 is any positive number such that S + 9 < lu A r/, p£ [1, 1/7) is chosen 
in order to fulfil the bounds 

1 p 

< z — — < r 



1 — 7 1 — 7p 

(p' is its conjugate exponent), and we utilized Proposition 12 .31 to go from (|4.19[) 
to (|4.20p . Notice that Proposition l2.3l applies. since the required constraint (12. 2p 
is satished. 

Thus, (|4~20|) shows that the map t m- e St e At Bu belongs to L r '(0,oo; [V(A* e )}'), 
r' being the conjugate exponent of r. From r G (1/(1 — 7), 00] we get r' G 
[1, 1/7), and (|4.18[) holds true with s = r' G [1, 1/7), as desired. 

Notice carefully that the range of the summability exponent s for the validity 
of (|4.18[) cannot be improved. In fact, owing to Proposition 12.31 the exponent r 
in the obtained estimate (|4.20j) is subject to the constraint (|2.2p , which implies 

p- (1 - TP) 
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while it is readily verified that 

P 1 

sup j- r = - , 

l<p<i. P~ (1 -IP) 7 

This confirms that (|4.18[) holds true if and only if s E [1, 1/7), thus concluding 
the proof. □ 

In order to pinpoint the regularity of $(/;)£?, we will employ the usual repre- 
sentation of the optimal state in terms of the initial state. It follows from (13. 5ft 
that 

<f>(-)x = (7 + LL*R*R)- x e A x , 

which becomes 

$(> = (I-LAT X L*R*R) e A 'x, (4.21) 

where A = I + L*R*RL is boundedly invertible on L 2 (0, 00; U) — an elementary 
consequence of the fact that A is coercive on L 2 (0, 00; U). 

Remark 4.8. We note that the representation of the inverse (I + LL* R* R)~ l 
which occurs in the formula (|4.21[) is easily derived by a direct (algebraic) com- 
putation. The inversion of an operator of the form / + SV in a Hilbert space 
setting is discussed in full detail in (30j Lemma 2A.1, p. 167]. 

If we take now x = Bu in (|4.21l) and formally rewrite the corresponding 
formula, we obtain the following representation for $>(t)Bu: 

e Apt Bu = <$>(t)Bu = e At Bu - [LA^ 1 L* R* Re A ' Bu] (t) , (4.22) 

where A -1 is required to make sense on the space L q (0, 00; U), or Li(0, 00; U) for 
some positive S, rather than on L 2 (0, 00; U). Indeed, given u 6 U, by Lemma l4"771 
we know that 

e A Bu £ LKO^oo; [D(A* e )}') . 
Then, owing to Proposition 12.91 the application of the operator L*R*R* yields 

[L* R* R*e A ' Bu] (t) e L q s {0, 00; U) , 

which holds for any S e [0,uj A 77 [. 

The question which then arises is the following. 

Question 4.9. Is the operator A = I + L*R*RL boundedly invertible on the 
function space L|(0, 00; U) ? 

It will become clear in the proof of Theorem 14.121 below that in contrast with 
previous theories A is not an isomorphism on L|(0, 00; U) (with a fixed 6), as 
we would expect. An intermediate useful result is the one given in the following 
Lemma. 

Lemma 4.10. The operator A = I + L* R* RL is an isomorphism in the space 
L 2 (0, 00; U), provided that S G (0, (J A 77) is sufficiently small. 



32 



Proof. We seek to solve uniquely the equation 

w + L* R* RLw = h, heLj(0,oc;U). (4.23) 

Let us denote by H the function space L^(0,oo;f7). Since H C L 2 (0, oo;U), 
there is a unique w G £ 2 (0, oo; £/) such that (|4.23l) holds. Multiplying (|4.23|) by 

e 5 *, we get 

e St w + e 5f L*R*RLw = e st h, 

which is equivalent to 

e st w + L* A _ 5 R*RL A+S (e st w) = e St h , (4.24) 

where we denoted by L A+ $ the input-to-state map for the control system y' — 
(A + S)y + Bv, y(0) = 0, namely 

t 



L A+S v(t):= / e^ A+5 ^ t -^Bv(s)ds: 



L* A _ S is defined accordingly. Thus, in order to simplify the notation, let us 
rewrite (I4.24[) as follows: 

e st w + L*_ 5 R*RL s (e St w) = e St h £ L 2 (0, oo; U) . (4.25) 
We now utilize the estimate 

\\L*_ 5 R*RL 5 - L*R*RL\\ ciL 2 (0 ^ u}) < c {r]A ^_ 5)2 , (4-26) 

which will be established in Lemma 14.111 below. The above implies that I + 
L*_ S R* RLs is also invertible in L 2 (0,oo;£7), with continuous inverse, for suffi- 
ciently small 5. Hence, the equation 

z + L*_ s R*RL s z = e st h 

has a unique solution z G L 2 (0, oo;U). Observe now that the function e~ St z 
belongs to H C L 2 (0, oo; U) and satisfies 

e" st z + L*R*RL(e- St z) = h . (4.27) 

Comparing Q4.23|) with the above gives, by uniqueness, e~ St z = w, so that 
w G H. This shows that (|4.23[) is uniquely solvable in H . The proof is completed 
once we establish the estimate (|4.26[) . This is accomplished in the following 
Lemma. □ 

Lemma 4.11. If 5 g]0,i)Aw[, the estimate (14.261) holds true. 

Proof. We have 

L*_ S R*RL S - L*R*RL = [L*_ s - L*]R*RL S + L*R*R[L 5 - L]. 

Consider now Lg — L. It holds for each u G L 2 (0, oo; U) 

[L s -L]y(t) = f e A{t ~ s) Bu(s)[e 5(t -^ - l]ds = 
Jo 

= f e s ^e A ^Bu(s)[l-e- s ^]ds, 
Jo 
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so that we can split 

[L s -L]y(t) = f e 5{t -^F(t-s)*u( S )[l-e- 8{t -^]d S + 
Jo 

+ f e s ^-^G{t-s)*u{s)[l-e- s ^-^]ds. 
Jo 



The first term can be estimated by 



/' 

Jo 



,6(t-a) 



F(t-s)*u(s)[l -e 



-5(t-s) 



< 



<cS [ (t-s 
Jo 



- aV—re-ii-W-s) 



\\u(s)\\uds, 



and it is straightforward to deduce that 



s) ]dt 



pOO pt 

/ / e 5 ^F(t-s)*u{s)[l- e - 
Jo Jo 

< cS [ a'-ie-^-^daWuh^.u) < C 
Jo 



dt 



< 



(77 - S) 2 -y 1 



lL 2 (0,oo;t/)- 



The second term is estimated as follows: fix ip S L 2 (0, 00; U) and set a = <5+T j A " , 

Q = SS+r^ = 5±a_ Then wg 

Jo"" (/ t ^ (t_S)G (*" s )* u ( s )[ 1 -^'"'l^.^)) dt = 



= J™ * u{s), e s{t - s ^ [1 - e-^-^Git - s)^(t)j dsdt = 

nt 
e-^-W-^Wu^Wu \\e^-^[l - e-^-^]\\G(t - s)^(t)\\u dsdt < 

poo 

/ \\e-^- s ^-M\L H o,uu)\\(t - -)e^G(t - OV'WIU^o,*^)^ < 
Jo 

ijc roc 

||e-^-*)( t -)«|| L2( o lt;to ||e a -G(-)^(t)||L=>(o,t^)*< 

pOO 

n , 5 J o \\e-^- s){t -M\L^ m \\m\Wdt. 



< cS 1 i 
'0 

26 

< c 



T) A ui — S J 
25 



< c 



From here it is a standard matter to deduce that 



< C 



{■n a w - £) 2 H U II L2 (°' 00 ^) 11^11^(0,00;!/) • 



This shows that 



/ e^ s) G(t - s)*u(s)[l - e-^-^ds 
^0 



< 



L 2 (0,oo ; y) 



< C 
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and summing up we obtain 



\\[L S - L]u\\ L 2 (0tOD . Y ) < C _ ^ 2 ||m||£2(o ;00;C 7) . (4.28) 



Next, consider L*_ s - L* . It holds for each y G L 2 (0, oo; Y) 

/oo 
B * e A>{T-t) [e - S[ -r-t) _ 1]y ( T)dT = 

X) f'OC 

F(r-t)[e- s( - T -V -i]y( T )d T + / G(j - t)[e- 5 ^ - l]y(r)dT 



From here, proceeding quite similarly to the preceding case, we get 

\\[L*- S - i*]j/||£=(0,oo;t/) < A ^_^ 2 blU 2 (0,oo;y) ■ (4.29) 

Finally, we can write 

[L*_ S R*RL S - L*R*RL)u = [L*_ s - L*]R*RL s u + L*R*R[L S - L]u 

and both terms can be easily estimated by (14.281) and (|4.29l) . The result follows. 

□ 

We now utilize Lemma T4. 101 to show that the operator A admits a bounded 
inverse A -1 which maps L|(0, oo; J7) onto L q s _ (0,oo;[/), for a suitable uq G 
(0,5). 

Theorem 4.12. The operator A = I + L*R*RL admits a bounded inverse 

A- 1 : L|(0, co; 17) — ► £^(0, oo; U) , 

im'i/l appropriate G (0,(5). 

Proof. We seek to solve uniquely the equation 

,g + L*R*RLg = h , (4.30) 

where /i G L q s (0, oo; {/), with arbitrary S G (0,uA 77). We follow an idea which 
has been employed in the study of the LQ-problem for parabolic-like dynamics; 
see, e.g., [301 Vol. I, Theorem 1.4.4.4, p. 40]. Since q < 2, using the action 
of both the operators L and L* , whose mapping increase the (time regularity) 
summability exponents, it is readily seen that there exists an integer no > 1 such 
that (L*R*RL) n °h G Lg(0, 00; U). Thus, we introduce the auxiliary equation 

v + L*R*RLv = (L*R*RL) n "h G L 2 (0, 00; U) . (4.31) 

Owing to Lemma T4. 101 possibly choosing 5 sufficiently small, (|4.31l) is uniquely 
solvable, yielding v G Lg(0, oo;U). Using once again that q < 2, it is easily 
verified that the obtained v belongs to Lg(0, 00; U) for any 9 < S. In fact, 

poo 

<»*\\v(t)\\ 9 udt= e-^- e ^\\e St v{t)\\ q udt (4.32) 
Jo 

( / e-M s - e ^-Wdt) (2 ~ 9)/ \ r\\e st v(t)\\ 2 udt) q/2 , 
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and there exists a constant C such that 

IMU*(0,oo;[/) < C \\ V \\l 2 6 (0,oo;U) ■ 

The function v resulting from (|4.3ip and the given h G W — for which the smaller 
9 still guarantees e e 'h G L q (0,oo;U) — will eventually produce the soughtafter 
solution g of equation (14.30[) . according to the following definition: 

no— l 

g= (-L*R*RL) j h + v. 

Indeed we have: 
(I + L*R*RL)g = 

"o- l n — 1 

= J2(- L * R * RL Y h + L*R*RL(-L*R*RL) j h + v + L*R*RLv = 

3=0 3=0 
no — 1 np — 1 

= j2(~ L * R * RL y h - Y(- L * R * RL y +lk+ (- L * R * RL ) n ° h= 

3=0 3=0 
no - 1 np 

= ^ {-L*R*RL) 3 h - Y^{-L*R*RL) i h + (- L* R* RL) n " h = 

3=0 i=l 

= h- (—L*R*RL) no h + (—L*R*RL) no h = h , 

which concludes the proof. □ 

We are finally able to show that the operator e Apt B substantially 'inherits' 
the regularity of e At B, except for a constraint on the exponent of the allowed 
exponential weights. 

Proposition 4.13 (Statement S6. of Theorem ll.5[) . For t > 0, the linear op- 
erator e Apt B is well defined as an operator from U into T>(A* e )\' and, in fact, 
provided 8 G (0, lj A rj) is sufficiently small, we have 

e s e Ap B G C(U,L p {0, oo;[V {A* ")]')) Vp G [1, — ) . (4.33) 

7 

Proof. Let u G U be given. We return to the representation (|4.22l) for e Apt Bu, 
and focus on its second summand. Starting from e At Bu, whose regularity is 
established in Lemma l4~7l we utilize (|2.14bl) of Proposition ^. 91 first, and invoke 
Theorem 14.121 next, thus obtaining — possibly choosing 8 < lj A rj sufficiently 
small — , 

hT 1 L* R* Re A ' Bu G L q e (0, oo; U) . 
In particular, there exists j3 < 9 such that 

hT 1 L* R* Re A ' Bu G L}(0, co; U) ; 

consequently, Proposition 12. 6f i) implies 

LA- 1 L*R*Re A Bu G Ll(0, oo; [V(A* e )}') Vr G [1, -) . (4.34) 
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Actually, the range of the summability exponent in (|4.34l) is larger; namely, 
the membership in (I4.34[) holds true for all r in a suitable maximal interval 
I D [1, 1/7), which is determined by the reciprocal relation between 1/(1 — 7) 
and q; see the statements (ii)-(iv) of Proposition !!?!)! Even an 'improved version' 
of the regularity in (|4.34j) , combined with the one in (14.18[) , yields anyhow 



e Ap Bu = e A Bu - [LA^L* R* Re A ' Bu}(-) £ L P (O, 00; [V(A* e )}') Vp £ [1 



for suitably small (3 £ (0, uj A 77), confirming (I4.33[) . □ 

The power of the (apparently weak) regularity result provided by Proposi- 
tion 14.131 is enlightened in the following Corollary. 

Corollary 4.14. The optimal state semigroup $(£) = e Apt is strongly differ- 
entiable on T>(A): more precisely, if x £ 13(A) the map t H> $(i)x = e Apt x is 
strongly differ entiable for almost any t > 0, with 

j + e Apt x = e Apt Ax - e Apt B B*Px , for x £ V(A) and a.e. t > 0, (4.35) 
and the equality holds true on [D(A* e )] r . In particular, the operator 

^-e Apt is continuous: V(A) — > L p s {0, 00; [V{A* e )}') Vp £ [1, -) (4.36) 
dt 7 

provided S £ (0,w A77AW1) is sufficiently small, and the following estimate holds 
true almost everywhere in (0, 00): 



d 

-e* 
dt 



A ~t^_„A P t 



< Mte-^WxWviA) 

(4.37) 

+ \\A- e e Apt B\ 



\c(u.y) \\ B * p \\c(v(A-),u) IMb(^) 



Proof. Let x £ T>(A). We preliminary note that e Apt x is strongly differentiable 
as an element of the dual space [X>(^4p)]'. In fact, if z £ V(A* P ), then z = A* P ~ 1 w 
with w £ Y and we may compute 

±(e Apt x,z) [v{Ap)V:V{Ap) = ^-(A p 1 e Apt x,w) Y = (e Apt x,w) Y 

= (e Apt x,A P z) Y = {A P e Apt x,z) [viA - p)] , MAp) = {x,e Apt Apz) Y 
= (x,Ape A " pt z) Y = {e Apt A P x,z) [v{Ap)YMAp) 
= (e Apt A[I - A~ l B B*P]x, z) [v{Ap)] , . v(Ap) , 

which shows that when x £ T>(A) 

3^-e Apt x = e Apt A[I - A~ l BB*P]x, for a.e. t > 0, (4.38) 

as an element of [D(A P )]' . In addition, since e Apt Ax £ Y, (14.38|) yields — still 
for x £ T>(A) and a.e. t > — , 

e A P t A [ A -1 B B *p} x = e A P t Ax _ e A P t A ^j _ A -l B B *p^ x g p(A P )]' . 
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Next, we observe that x G T> (A) implies x G T>{A e ) for all e G (0, 1), which 
ensures B*Px G U in view of Theorem GOD Then BB*Px G [D(A*)}' and 
also A A^BiTPx G [£>(A*)]', because AA -1 coincides with the identity oper- 
ator on [©(A*)]'. We now recall (|4.33[) from Proposition 14. 131 which establishes 
e Apt B B* Px[D(A* e )]' , along with the regularity (in time) result 

e s e Ap BB*Px G L p (0, oo; [X>(A* £ )]') Vp G [1, — ) , (4.39) 

7 

valid for sufficiently small 6 > 0. On the other hand, since the semigroup e Apt 
is exponentially stable (with the estimate (jl.lll) . we find 

e s 'e Ap ' Ax G L s (0, oo; Y) VsG[l,oo], (4.40) 

for any S < u)\. In view of the memberships (|4.39l) and (|4.40l) we see that 

e Ap ' A{I—A~ 1 B B*P]x = e Ap Ax~e Ap B B*Px G L p s (0, oo; [V{A* e )]') , (4.41) 

provided 6 is sufficiently small. 

Thus, (14.41[) shows that the derivative in (|4.38[) — a priori taking values on 
[D (A* P )] ' — coincides with e Apt Ax - e Apt A- 1 BB*Px G [V{A* e )]' for a.c. t > 
(We recall that the inclusion [V(A* P )]' D [D(A* e )]'— which holds true provided 
e < 1 — 7 — is the dual statement of (|4.8[) of Proposition 14.51 whereas here e can 
be taken arbitrarily small.) Therefore, (|4.35p actually makes sense on [©(A* 6 )]' 
for a.e. t > 0, and the validity of (|4.36[) is established, provided S is sufficiently 
small. 

Finally, natural estimates for each summand in the right hand side of (|4.35|) 
produce the bound in (|4.37[) . thus completing the proof. □ 

4.3 Well-posedness of the ARE 

We begin with the statement of a Lemma which provides boundedness of the op- 
erators A*P and A* P P on appropriate spaces (T> (A) and V(Ap), respectively), 
properties which constitute a prerequisite for well-posedness of the ARE. Al- 
though the proof is fairly standard, it is given below for the reader's convenience. 

Lemma 4.15. The following statements pertain to the optimal cost operator 
P. 

(i) A*P G C{V{A P ),Y), with 

A*Px = -R*Rx - PA P x Vie G V(A p ) ; (4.42) 

(ii) A* P P G C(V(A),Y), with 

ApPx = -R*Rx - PAx Vx G V{A) . (4.43) 

Proof, (i) Let x G V(Ap). Write the formula (|1.12j) which defines the Riccati 
operator, that is 

/•oo 

Px= I e A ' t R*R<I>(t)xdt, 
Jo 
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and integrate by parts in £, thus obtaining 

Px= A*~ 1 A*e A ' t R*R<S>(t)xdt 
Jo 



= A*- 1 e A * t R*R$(t)x - I A*~ 1 e A ~ t R*R<f>(t)A P xdt 



t — OO 



t=0 

oo 



I * — 1 ^A* t T)* 







p CO 

= —A*~ x R*Rx - A*^ 1 / e A " t R*R<£(t)A P xdt 
Jo 

= -A*~ 1 R*Rx - A*~ 1 PA P x . 

The above identity shows that P maps V(A P ) into T>(A*), and also that (|4.42p 
holds actually in Y, since R*Rx — PApx G Y for any x G 2?(^4p); the bound- 
edness of A*P immediately follows. 

(ii) We write 

(Px,z) Y =( [ e A * t R*Re Apt xdt,z\ =[ (e A ' t R*Re Apt x,z) Y dt 
^Jo ' Y Jo 

= [ (x,e Ap ' t R*Re At z) Y dt= (x, f e Ap ' t R*Re At zdt) =(x,P*z) Y . 

Jo V J 1 Y 

Since we know that P — P* , we deduce the alternative formula 

/■oo 

Px= e A * pt R*Re At xdt, x EY . (AAA) 



Thus, the proof of (I4.43[) follows almost precisely as in the proof of (|4.42l) . 
bringing about the equality 

Px = -A* P ~ 1 R*Rx - A* P ~ x PAx xeV(A), 

which confirms that A* P P is a bounded operator on T>(A) with values in Y, and 
hence (|4.43l) holds in Y. 

□ 

We are finally ready to show that the optimal cost operator is a solution of 
the the Algebraic Riccati equation corresponding to Problem II .31 

Theorem 4.16 (Statement S7. of Theorem 11.51) . The optimal cost operator P 
defined in (|1.12[) satisfies the following regularity properties: 

P e C{V(A P ), V(A*)) n £(V(A),T>(A* P ) . 

Moreover, P is a solution to the Algebraic Riccati equation 

{Px, Az) Y + (Ax, Pz) Y - (B*Px, B*Pz) v + (Rx, Rz) z = 

for any x, z G T> (A), 

which reads as 

(A*Px,z) Y + (x,A*Pz) Y - (B*Px,B*Pz)u + (Rx,Rz) z = 

when x, z G T>(Ap). 



(4.45) 



(4.46) 
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Proof. The proof splits into two parts. First, we establish the validity of (I4.45[) 
for x, z £ D(A). Beside Lemma [4.151 Corollary 14.141 will provide the crucial 
tool. Next, we show the validity of (|4.45l) with x, z £ T>{Ap). We will use once 
again Lemma 14.151 along with the intrinsic representation (14.31) of Ap in terms 
of the operator T. 

1. We recall the alternative representation (I4.44[) of the Riccati operator 
P obtained in the previous Lemma, and for x, z £ Y write the inner product 
{Px, z)y as a function of t: 

j* oo />oo 

(Px,z)y= {Re Apt x,Re At z) z dt = {Re Api - T ~ t] x, Re A{T - t] z) z dr. 

(4.47) 

Taking now x, z £ T>(A), in view of (I4.35|) of Corollary 14. 141 we differentiate 
both sides of the obtained formula with respect to t to find 

= — {Px,z) Y = -{Rx,Rz) z 

/CO 
(e A ^A P x,R*Re A ^z) mMt)yMA , t) dr (4.48) 

/OO 
(Re A ^ T ' t ^x,Re A ^ T - t '>Az) z dT 

/OO 
^p{r-t)Ax, R*Re A( - T ~ t h) Y dT 

/OO 
(e A ^BB*Px,R*Re A ^z) mA , nyv{A , n dr 

/OO 
(Re A ^-^x,Re A ^Az) z dr 

= -(Rx, Rz) z - (Ax, Pz) Y + (B*Px, B*Pz) v - {Px, Az) Y , 

so that 

(Ax,Pz) Y + (Px,Az) Y - {B*Px,B*Pz)u + (Rx,Rz) z =0, x, z £ V(A) . 

(4.49) 

Note carefully that the integrand in (|4.48[) makes sense a priori as a duality 
pairing on T>(A* e ) (which is summable because of (|4.36p ). whereas all the sum- 
mands in (|4.49p actually make sense as inner products in the spaces Y, U or Z, 
since in view of Theorem 13.51 B* P is bounded (a fortiori) on T>(A). Therefore, 
P solves the Algebraic Riccati Equation (|4.451> on V(A). 

2. We preliminarly recall that if a; £ V(Ap), then Apx = A (x + Fx), 
where the operator T — which is defined by the equivalent weak limits (|4.6[) and 
(|4~7)) — coincides on V(A P ) with -A~ 1 BB*P; see Remark Ol 

Given now x, z £ T>(Ap), we compute 

(A*Px, z) Y + (x,A*Pz)y = (by Lemma |4TT5|) 

= -(Rx,Rz) Y - (PA P x,z) Y + (x,A*Pz) Y = 

= -{Rx, Rz) Y - (PA{x + Tx),z) Y + (x, A*Pz) Y = 

= -{Rx, Rz) Y - {x + Fx, A*Pz) y + (x, A*Pz) Y = 

= -{Rx,Rz) Y - {Tx,A*Pz)u . (4.50) 
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On the other hand, we readily have 

(Tar, A*Pz) Y = -(A^B B*Px, A*Pz) Y 

= -{B B*Px, P«)p> (A . )]v > ( A.) = -(B*Px, B*Pz) Y , 

with B*P G C(V(A P ),U). Thus, inserting (|43Tj) in ([430)1 we find that P 
satisfies the Algebraic Riccati equation (|4.46[) for any x, z G 2?(Ap), with _B*P G 
£(D(ylp), [/), thus concluding the proof. □ 



5 Illustrations 

In this section we give a significant illustration of the applicability of the infinite 
time horizon optimal control theory provided by Theorem 11.51 The boundary 
control problem under examination is the thermoelastic plate model studied in 
[To] and [Tj. Recall that this specific PDE problem constituted the prime mo- 
tivation for the introduction of the novel class of control systems characterized 
by the abstract assumptions listed in [2j Hypotheses 2.2]. 



5.1 A thermoelastic plate model with boundary thermal 
control 

We consider a classical (linear) PDE model for the determination of displace- 
ments and the temperature distribution in a thin plate; see [521 121] • Let Jl be a 
bounded domain of M 2 , with smooth boundary T. The PDE system comprises a 
Kirchhoff elastic equation for the vertical displacement w(x, t) of the plate and 
the heat equation for the temperature distribution 9(x,t). The plate equation 
is supplemented with clamped boundary conditions, whereas a control action 
on the temperature, represented by the function u(x, f), is exercised through T. 
Thus, the PDE problem reads as follows (the constant p is positive, v denotes 
the unit outward normal to the curve T): 



w tt - pAwtt + A 2 w + A6 = 
e t -A9- Aw t = 
^ = 1^ = 

e = u 

w{Q,-) = w°, w t {0,-) = w 1 ; 6(0,-) = 



in SI x (0, oo ) 
in Q x (0, oo) 
on T x (0, oo) 
on T x (0, oo) 



(5.1) 



CI. 



With (|5.1I) we associate the natural quadratic functional 



(\Aw(x, t)\ 2 + \Vw t (x, t)\ 2 + \9(x, t)\ 2 ) dx dt- 



\u(x,t)\ 2 dsdt (5.2) 



to be minimized overall u G L 2 (T x (0, oo)), where w solves the boundary control 
problem (|5.1[) . Note that when u = the functional (|5.2[) is nothing but the 
integral over (0, oo) of the physical energy E(t) of the system. 

We recall that the finite time horizon optimal control problem for the con- 
trolled PDE system (|5.1[> was first studied in [16] , and then fully solved according 
to the novel abstract theory set forth in [5] . The preliminary PDE analysis car- 
ried out in [16] . combined with the key boundary regularity result established 
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in [TJ Theorem 1.1.] provided the proof. 

Here we aim to complete the study of the associated optimal control problems 
including the infinite time horizon case. Specifically, we will show that the 
model under investigation fits as well into the abstract framework designed by 
Assumptions 11.11 and 11.41 thereby ensuring the applicability of Theorem 11.51 

We already know that the boundary control problem (|5.ip can be recast as 
an abstract control system of the form (|1.1[) in the state variable y — (w, Wt,0), 
with appropriate (and explicitly derived) dynamics and control operators (A, B); 
see [IHl § 2] for all details. A thorough analysis of the semigroup formulation of 
the uncontrolled problem is provided by [27] , where the predominant hyperbolic 
character of the coupled PDE system (in the case p > 0) was first pointed out. 
The state and control spaces are given by 

Y = Hg(Q) x i?o(ft) x L 2 {Vt) , U = L 2 (T), 

respectively. The plan is thus to check the complex of requirements contained 
in Assumptions I1.1H1.4I 

Verification of Assumption [TTTJ The validity of the basic Assumption ll.il 
has been already discussed in [T5], yielding the explicit statements of [16, Propo- 
sition 2.1.]. We recall from [161 Remark 2.2] that well-posedness of the uncon- 
trolled model was proved in |27j , while an easy computation provides bounded- 
ness of the linear operator A~ 1 B; see (2.23) in [T|5]. Instead, the exponential 
stability of the underlying semigroup — by far a more challenging issue — was 
established in [S]. 

Verification of Assumptions ITT41 We must verify that all the requirements 
listed in Assumption 11.41 are fulfilled. Accordingly, we recall from [TJ § 5] that 
given z — (iu ,^ 1 ,^ ) e V(A*), one has 



A-t 99 
B e *z - j- 
ov 



(5.3) 



where now 8(x, t) is the thermal component of the solution y(t) = (w(t), Wt(t), 0(t)) 
to an initial/boundary value problem which comprises the same thermoelastic 
system of (|5.ip . yet with homogeneous boundary conditions and with a slightly 
different initial condition, that is y(Q) = (w° , —w 1 ,8°). The change of sign 
on a component of initial data is not influential, and it justifies the estimates 
performed on the solution to the original PDE problem (|5.1[) with u = 0. 

For the reader's convenience we record from [TJ § 5] the essential steps of 
the computations leading to the sought decomposition (|1.5D of B*e A t yo, with 
suitable F and G that will be shown to satisfy the series of assumptions in (|1.4j) . 
1. Notation. In the formulas below the symbol — Ad denotes the (un- 
bounded) linear operator which is the realization of the Laplace operator A 
in H = L 2 (fl), when supplemented with (homogeneous) Dirichlet boundary 
condition, i.e. 

A D w = -Aw, w e V(A D ) = H 2 {Q) n H£(Q) . 

It is well known that —Ad is the generator of a strongly continuous semigroup 
e — At>t j n jj which moreover is analytic. The fractional powers A D are well 
defined for all a S (0, 1), and there exist a positive constant lol and constants 
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L a > 1 such that the following estimates hold true: 



e 



-ui L t 



W D e- A ^\\ c{H) <L a — aG[0,l]. (5.4) 

The related (positive) operator M. = I + pAo is employed in the abstract 
formulation of the elastic equation and hence will occur in the computations 
below. Instead, A will denote the realization of the bilaplacian A 2 in O with 
homogeneous clamped boundary conditions. 

Finally, let D be the map which associates to any function in L 2 (T) its 
harmonic extension in il. Classical trace theory ([33]) yields 

D continuous : L 2 (T) — ► H 1/2 (n) C V(A^ a ) , < a < i , 
which implies 

A^ a D continuous : L 2 {T) — > L 2 (il) , < a < - . (5.5) 

In addition, the following well known result will be used throughout: 
dh 



D*A D h= — 
av 



r 



fe 6 ff3/a+'(n)nflj(n), <t>o ; (5.6) 



see, e.g., [3D1 Lemma 3.1.1, p. 181]. 

2. Explicit decomposition. We take the explicit expression of 0{t) and 
utilize (15.61) to rewrite 



^ =D*A D 9{t) = D*A D e- ADt 6°-A D / e-^-'WaJda 
av L ,/ n 



we next integrate by parts, obtaining first (since according to the clamped 
boundary conditions =0on T): 

m_ rt 

dv 



= D*A D e- ADt 9° + D*A D e- ADt w 1 -D*A D [ e- AD{t - s *>w tt (s) ds . 



Fi(t)j/ F a (t)vo 

V>(*;yo) 

Thus, we utilize the elastic equation w u — — M~ 1 Aw(s) + M^ 1 AdO{s) to find 
V>(t;yo) = -D*A D I e- AD{t - s) M- 1 Aw{s)ds 



V>i(*;yo) 



D*A D f e^ AD ^ s ^M~ 1 A D e{s)ds 
Jo 



= D*A D f e - AD( - t - s) M- 1 A D Aw(s)ds 
Jo 



mt)y ° (5.7) 



■0ii(*;?/o 
ft 



—D*Ajj / e-^-^M-'Ao Aw(s)\ r ds -F 3 (t;y ) 
Jo 



i>i2(t;yo) 
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where the further splitting of ifii(t; yo) = i/Ju(t; yo) + ipi^it; yo) is a consequence 
of 

M~ x Aw = -M~ 1 A D (Aw - DAw\ r ) ; 

see [ini §5, formula (5.9)]. We will see that the integral ipn(t',yo) eventually 
contribute to the term F(t)yo, while it was shown in [16] that ^12(^2/0) does 
satisfy Assumption ll.4f ii) and (iii)(a); hence, V'i2(i; yo) is identified with G(t)yo. 
Summarizing, we found the following decomposition: 

89 4 

— =X> i (i)j/ + G(i)yo, 
1=1 

where 

Fi(t)y = D*A D e- A ° t 6° , F 2 (t)y = D*A D e- A ° t w 1 

F 3 (t)y = -D*A D f e- AD ^M- 1 A D B{s)ds, 
Jo 

F A {t)y Q =D*A D [ e- AD ^ t -^M- 1 A D Aw(s)ds , 
Jo 

and 

G(t)y = D*A D [ e- AD( - t -^M- 1 A D Aw{s)\ r ds. (5.8) 
Jo 

3. Estimates. We proceed to check the validity of the first among Assump- 
tions [T31on either term Fi(t)y$, i = 1,...,4. Employing the various properties 
recorded above — in particular, using repeatedly the analytic estimates (|5.4I) — 
and that 

\\(w(t),w t (t),9(t))\\Y = \\e At (w ,-w\d°)\\ Y <Me-^\\y4 Y , 
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we find 

\mt)y Q \\ <\\D*A¥-°\\L 3/ t +(TW ^\\e \\ <C 1)ff r 3 / 4 - CT e— 1y Q ||; (5.9) 
\\F 2 (t)y \\ = \\D*A^ 2 \\ We-^A^W 

< C 2 , CT 4 max{t- 3 / 4 - ff , l}||y || ; (5.10) 

< |p*^ 4 - CT ]|i 3 /4 + . \\M- X A D \\ J* *~" L *~i +a Me— ||yo|| d S 
<C3, CT t 1/4 - ff e^ 4 bo ||, cJ = min{ W , Wi }; (5.11) 

||F 4 (i)y || < \\D*A^-°\\L Vi+a j* ^wT+a A D \\\\ W {s)h, u ds 

||Z?*< 4 -l£ 3 /4 +CT jT* ll-M-M^H Me— ||y || & 

<C4, CT t 1/4 ^e^ 4 1|2/ 1|, ZJ = mm{uj,uj L }. (5.12) 

Thus, if a is fixed, the bounds obtained in (|5.11l) and (15 . 12[) imply that 

3C, 77 > : ||^(t)2/o|| <Ce-" 4 ||2/ ||, z = 1,2. (5.13) 

(We note that, more precisely, for any r\ < min{o;,u;L} there exists a constant 
C v > such that the estimate in (|5.13[) holds. However, the formulation (|5.13|) 
will suffice.) 

Consequently, with a fixed a the estimates (|5.9p and (I5.10[) combined with 
(|5.13p show that Assumption ! 1 .4f i) is satisfied with 7 > 3/4 and rj < min{oj, u>l}- 

Finally, that the component G(t)yo in (|5.8p fulfils the requirements of As- 
sumptions |1.4f ii) and (iii) (a) was proved in |16| . The more challenging As- 
sumption [13{hi) — whose PDE counterpart is a suitable regularity result for the 
boundary traces of 9t — has been established in [TJ. 

Therefore, provided the observation operator R meets Assumption !! ,4f iii) (b) . 
all the hypotheses of Theorem 1 1.5 1 are satisfied; consequently, all its statements 
follow. Indeed, in the present case the observation operator which occurs in 
the functional (|5.2[) is the identity, and it is not difficult to ascertain that As- 
sumption [HlKiii)(b) holds true if e is taken sufficiently small; see [TH1 §3] and [H 
Remark 2.5]. 

A Appendix 

Proof of Proposition [X7| We shall use the representation (14.21)) for 3>(t)x. Let 
x G T>(A e ), with e given by Assumption 1 1 ,4f iii) . By standard semigroup theory 
we know that for any 5 < lu, 

e A ie^(0,oo;D(A ( )) Vpe[l,oo]; (A.l) 
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then, owing to Assumption 1 1 .4f iii) (b) . we obtain as well 

R*Re A x G L p s (0,oo;T>(A* e )) Vp G [1, oo] . 

In particular, R*Re A 'x G L^(0, oo; V(A* € )) with p > 1/(1-7) an d Proposi- 
tion [2?7]lv) establishes 

L*R*Re A x G L|°(0,oo;f7) . 

We now make the following Claim. 
Claim A.l. The operator A = I + L*R* RL admits a bounded inverse 

A- 1 : Lf(0, 00; U) — > L^ ff0 (0, 00; U) , 

wz</i appropriate <7o G (0,5). 

Assuming that Claim |A~T1 is valid, it follows that 

A~ 1 L*R*Re A 'x G Lg°(0, 00; £/) for some 9 < 5 < lo A rj. 

Thus, using now Proposition 12. 6( v) we see that 

LA- 1 L*R*Re A x e Lg°(0, 00; £>(A e )) c Ljg(0, 00; £>(A e )) (A.2) 

for any (3 such that < j3 < 9 and for all p G [2, 00]. 

The regularity in (|A.2[) , combined with the one in (| A. 1|) . shows that p.22p 
holds true. Then, using the boundedness of the gain operator B*P on V(A e ) in 
the feedback representation p. 241) of the optimal control, we easily obtain that 
(|3.22p implies Q3.23p . thus completing the proof. □ 

Proof of Claim \A.1[ We proceed pretty much in the same way as in the proof of 
Theorem l4.12l The tools which play a major role are, once more, the smoothing 
properties of the operators L and L*, as well as the inclusions L r a (0, 00; U) C 
L%(0, 00; U) for (3 G (0, a) and r < s (including s — +00). 

1. We seek to solve uniquely the equation 

g + L*R*RLg = heL^(0,oo;U), (A.3) 

with arbitrary S G (0, u A rf). An elementary calculation shows that h G 
Lg(0, 00; U) as well, for any < 9 < 5. Thus, since by Lemma 14.101 A is 
boundedly invertible in Lg(0, 00; U) provided 9 is taken sufficiently small, there 
exists a unique function g G Lg(0,oo;{7) such that (|A.3|) is satisfied. We will 
show that — possibly taking a smaller 9 — in fact 

5 eL e °°(0,oo;t/). (A.4) 

2. Preliminarly, we prove that 

3ni G N : {L*R*RL) ni continuous: Lj(0, oo; U) -> Lf{Q, 00; U) . (A.5) 

Let / G Lg(0, 00; U). Since the successive application of the operators L and L* 
improve the (time regularity) summability exponents, there exists an integer n 
such that 

{L*R*RLff G L q e '(0,oo;U), 
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which in view of Proposition^^) implies e 9 ' L(L*R* RL) n f e C b ([0, oo); V(A e )). 
Then, e e R*RL{L*R*RL) 1T f £ C 6 ([0, oo); X>(v4* £ )) and by Proposition [ST^v) 

{L*R*RLf +1 f = L* R* RL{L* R* RLf f € L^°(0, oo; CT) . 

Thus, (|A.5I) is satisfied with n\ = n + 1. 

3. We return to (|A.3[) and argue as follows. The solution g of (|A.3|) is such that 

g = h — L*R*RLg ; (A.6) 

then, if m = 1 in (IA.5I) . since /i and L*R*RLg in (jA.6[) both belong to 
Lg°(0, oo; [/), we immediately obtain (|A.4|) . The case n\ > 1 is treated as 
follows. 

We apply the operator L*R*RL* to both members of (|A.6|) . thus obtaining 

L*R*RL*g = L*R*RL*h - (L*R*RL*) 2 g . (A.7) 

Observe now that L*R*RL*h 6 L^(Q, oo; t7), as the given regularity of h is 
maintained by the application of L*R*RL; if, in addition, rii = 2 in ([A"Jj]> . the 
identity (|A.7p yields L*R*RL*g £ £g°(0, oo; [/). Using the obtained regularity 
for L*R*RL*g in ([P]) . we find that 5 G Lg°(0, oo; U), that is (jX4|) . 

If, instead, n\ > 2, inserting (|A.7[) into (|A.6[) we find the novel identity 

g = h - L*R*RL*h + (L* R* RL*) 2 g . (A.8) 

To pinpoint the regularity of the summand (L*R*RL*) 2 g, we return to (|A.7|) 
and apply the operator L*R*RL* to both members. Plugging the obtained 
expression for (L*R*RL*) 2 g into (1A.8[) yields 

2 

3=0 

The iteration of this argument eventually yields 

n i — l 

9=^2 (-^y(L*R*RL*yh+(L*R*RL*) ni g, 

3=0 

where both summands on the right hand side belong to Lg°(0, oo; U), which 
establishes (IA.4I) and concludes the proof. 

□ 
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